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Abstract We study products of arbitrary random real 2x2 matrices that 
are close to the identity matrix. Using the Iwasawa decomposition of SL(2, R), 
we identify a continuum regime where the mean values and the covariances 
of the three Iwasawa parameters are simultaneously small. In this regime, the 
Lyapunov exponent of the product is shown to assume a scaling form. In the 
general case, the corresponding scaling function is expressed in terms of Gauss' 
hypergeometric function. A number of particular cases are also considered, 
where the scaling function of the Lyapunov exponent involves other special 
functions (Airy, Bessel, Whittaker, elliptic). The general solution thus obtained 
allows us, among other things, to recover in a unified framework many results 
known previously from exactly solvable models of one-dimensional disordered 
systems. 
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1 Introduction 

Products of random matrices arise quite naturally in the context of one- 
dimensional or quasi-one-dimensional disordered systems [1-4] . Such products 
typically grow exponentially with the number of factors. The associated growth 
rate 7 is called the Lyapunov exponent of the product of random matrices. 

This paper deals with products of 2 x 2 random matrices with real entries 
and unit determinant, i.e., elements of SL(2,R). The relevance of products of 
2x2 transfer matrices to one-dimensional disordered systems was already un- 
derlined in the pioneering works by Dyson [5] and Schmidt [6] . In the context 
of quantum mechanics in a one-dimensional random potential, the exponen- 
tial growth of random matrix products is at the basis of the phenomenon of 
Anderson localization. The Lyapunov exponent is interpreted as the inverse 
of the localization length [7]. In classical statistical mechanical models, such 
as spin chains with random couplings and/or random fields, the Lyapunov 
exponent is the reduced free energy per site [2,3]. 

Besides perturbative weak-disorder expansions or numerical investigations, 
and the Lloyd model, with its Cauchy distribution of random site energies, 
which is solvable in any dimension [8], the cases where the Lyapunov expo- 
nent has been obtained in analytical form are rather few, even in the simplest 
situation of 2 x 2 matrices. These exact solutions concern classical disordered 
harmonic chains [5,9], variants of the disordered Kronig-Penney model [10, 
11], the tight-binding Anderson model [12], and classical and quantum spin 
chains [13-17]. See also the monographs [1,3,18,19] and, for more recent ex- 
amples, [20-23]. The success of the calculation in all those cases relies upon 
special properties of the distribution of the random matrices, with a single 
random variable whose distribution is simple (i.e., exponential or power-law). 

The goal of the present work is to gain greater insight into the generic 
case by considering distributions concentrated around the identity matrix. 
Our approach differs from perturbative approaches such as those of Derrida 
et at [24] and Sadel and Schulz-Baldes [25], in that we solve exactly a limiting 
form of the integral equation for the Lyapunov exponent. As we shall see, 
by using the Iwasawa decomposition of SL(2,R) into compact, Abclian and 
nilpotcnt subgroups, one can define a continuum regime in which the Lyapunov 
exponent exhibits a scaling form. Our main finding is that this scaling form 
can always be expressed in terms of the logarithmic derivative of a special 
function (Airy, Bessel, Whittaker, elliptic, hypergeometric) . We thus recover 
in a unified framework many results known previously from exactly solvable 
models of one-dimensional disordered systems, and obtain several novel results. 

1.1 Contents of the paper 

The paper may be divided into four parts: 

1. The first part is the remainder of this introduction, in which we give an 
overview of our results and techniques. After some basic definitions and 
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some heuristics, we identify the continuum regime of interest, and present 
the key equations; these are Equations (1.32) and (1.42) for a certain in- 
variant density and its Hilbert transform in the continuum regime, and 
Formula (1.39) for the Lyapunov exponent. The introduction ends with 
a concrete calculation that illustrates our general approach in a simple 
particular case. 

2. The second part, Sec. 2, contains a careful derivation of the key equations. 

3. The third part is the core of the paper. In Sec. 3, we show how the character- 
istic exponent — essentially a complexification of the Lyapunov exponent — 
can be systematically expanded in powers of the covariances of the Iwa- 
sawa parameters. We also compute the first two non-trivial terms of this 
weak-disorder expansion. The following sections are concerned with the 
exact calculation of the characteristic exponent in several situations, in 
order of increasing difficulty. Sec. 4 is devoted to the cases where only 
one of the Iwasawa parameters is random; for physical reasons that will 
emerge, the compact, Abelian and nilpotent cases that arise are referred 
to as distance, supersymmetric, and scalar disorder respectively. The case 
where there is both scalar and supersymmetric disorder, but no distance 
disorder, is treated in Sec. 5. We consider in Sec. 6 the situation where 
the three parameters are uncorrelated and have zero mean. Sec. 7 deals 
with the completely general case, where the characteristic exponent may 
be expressed in terms of Gauss's hypergeometric function. 

4. The fourth part brings out the relationships between our results and vari- 
ous aspects of the theory of disordered systems. Sec. 8 develops the connec- 
tion with Brownian motion in the Poincare hyperbolic half-plane. Sec. 9 
discusses the interpretation of our results in terms of a disordered quantum- 
mechanical model. Finally, our findings are briefly summarised in Sec. 10, 
while an appendix is devoted to the derivation of Equation (1.42). 



1.2 The Furstcnbcrg formula 

We address the problem of computing the Lyapunov exponent 

:= lim M (L1) 

n— >oo ft 

of the product 

n n :=M n M n _ 1 ---M 1 , (1.2) 

where the M n are independent random matrices in SL(2,R) with the same 
probability measure m. Here, E(-) stands for the expectation with respect 
to m, and the notation | • | refers to the Euclidean norm in R 2 and to the norm 
on 2 x 2 matrices induced by it. 
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Our starting point the Furstenberg formula [1,26,27]: 

M 



7 



= r p(dz) i 

J-oo JS 



m(dM) In 



SL(2,R) 



(1.3) 



In this expression, z is the projective coordinate, i.e., the reciprocal of the 
slope associated with a direction in M 2 , and p is the probability measure on 
the projective line (the set of all directions) which is invariant under the action 
of matrices drawn from m. We emphasize that, whereas m may be considered 
as given, one essential difficulty of the calculation is that p must be found. In 
the particular case where p has a density, i.e., 



(1.4) 



it may be shown that the unknown density obeys the integral Furstenberg 
equation — often referred to as the "Dyson-Schmidt equation" in the physical 
literature: 



L 



m 



SL(2,R) 



{AM) (foM-^iz) 



dM^jz) 
dz : 



(1.5) 



where M. 1 (z) is the inverse of the Moebius (linear fractional) transformation 

m\\z + mi2 



M(z) := 



m 2 iz + m 2 2 



associated with the matrix 



M := 



mil 
m 2 i 



mi2 
m 2 2 



(1.6) 



(1.7) 



As pointed out earlier, there is no systematic method for solving (1.5); exact 
solutions are limited to very specific forms of the distribution m of the random 
matrix M. 



1.3 A one-dimensional disordered system 

Our approach builds to a large extent on the intimate connections between our 
topic and the theory of one-dimensional disordered systems. In this subsection, 
we describe a simple quantum-mechanical model that brings out the most 
relevant of these connections. 

Consider the Schrodinger equation on the positive half-line 

-tp"(x) + V(x) i){x) = E ij){x), (1.8) 

where the potential V consists of arbitrary point scatterers located at the 
points 

< X! < x 2 < ■ ■ ■ (1.9) 
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This means that the potential vanishes everywhere, except at the positions x n 
of the scatterers. The action of each scatterer on the wave function ip is encoded 
in the linear rule 



, , - i = B n I Y .}^> (1.10) 

which relates the values of ip(x) and of its derivative ^'(x) to the left and 
to the right of the scatterer. The requirement that the resulting Schrodinger 
operator have a self-adjoint extension translates into the condition [28,29] 

e i9 B n g SL(2,R) (1.11) 

for some Set, and there is no appreciable loss of generality in taking 9 = 0. 
We then have 

where the matrices in the product are given by 

v° 7f/v sina « cosa„ y "^o Vfcy v y 

and belong to SL(2,M). Here, k = \f~E is the momentum of the particle, 
a n = k£ n , and £ n = x n+ i — x n is the distance between neighbouring scatterers. 
The most familiar example is that of the standard (scalar) delta-scatterer 
introduced in the deterministic case by Kronig and Penney [30], and studied 
in the random case by Frisch and Lloyd and others [31,32]: 

Bn=(l (1-14) 

Frisch and Lloyd investigated the density of states of this model by making 
use of the Riccati variable 

*> - tr ('■«) 

viewed as representing the "position" at "time" a; of a particle moving along 
the real axis. They observed that the integrated density of states per unit 
length is the (negative of the) stationary current of particles or, what is the 
same, the reciprocal of the mean time that the Riccati variable takes to make 
a complete journey along the real axis. 

In his later study of the Frisch-Lloyd model, Halperin [33] made two further 
contributions. Firstly, he considered the particular limit of the model where 
the scatterers become infinitely dense whilst their strength becomes infinitely 
weak; we call this the continuum regime. Secondly, Halperin recognised that 
it was not strictly necessary to know the stationary distribution of the Riccati 
variable in order to find the integrated density of states; he showed how it could 
be obtained indirectly by considering a certain transform of the stationary 
density. While these early works were concerned with the integrated density 
of states, it was eventually recognised that the Lyapunov exponent and the 
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integrated density of states are tightly related. They can indeed be viewed 
as the real and imaginary parts of a function analytic in the complex energy 
plane — a property that results in the Herbert- Jones-Thouless formula [34,35]. 

The approach and the results we present in this paper can be thought of 
as extensions of Halperin's ideas to the case where the delta-scatterer (1.14) 
is replaced by the more general point-scatterer 

B »=(T X)(l T)- <>•»> 

with a scalar component of intensity u n and a supersymmctric component of 
intensity w n . Comtet et al. called this point-scatterer the double impurity [21]. 
The relevance of our results to this quantum-mechanical model will be dis- 
cussed in Sec. 9. 



1.4 The complex characteristic exponent 

Our first task is to identify the quantity that is, for the general product of 
matrices 77„, the counterpart of the integrated density of states (or current) 
in the disordered quantum model of the previous subsection. To this end, we 
integrate (1.5) with respect to z: 

Ely f(t)dt\=j. (1.17) 

The constant of integration j is completely determined by the requirement 
that / be an invariant probability density. 

We can gain some insight into the significance of the quantity j by intro- 
ducing a "Riccati process" {z n } defined by the random recurrence relation 

z n+1 = M n+ i(z n ). (1.18) 

Denote by /„ the probability density of the random variable z n . It then follows 
easily from this recurrence relation that 

f n+1 (z) = E { ^ M 2 [Z) % ° M- 1 ] (*)) (1.19) 
and, after subtracting f n {z) from both sides, 

fn+l(z) - fn(z) = E [/„ o M- 1 ] (z) - f n (z)j . (1.20) 

This discrete-time Fokker-Planck equation expresses the fact that the itera- 
tion (1.18) leads to a redistribution of the probability density of the projective 
variable. By introducing the quantity 

Jn (z):=E[ / f n (t)dt , (1.21) 
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we can write (1.20) in the "conservation form" 

fn+l(z) ~ fn(z) = -^-jn(z). (1.22) 

Therefore j n (z) has an obvious interpretation as (the negative of) the "prob- 
ability current" induced by iteration. Under mild assumptions on the distri- 
bution m of the matrices in the product, Furstenberg's theory asserts the 
existence of a stationary distribution of the projective process, and of a corre- 
sponding steady current: 



f(z) = lim f n (z) and j = lim j n (z). 



(1.23) 



This motivates the following choice for the definition of the complex char- 
acteristic exponent associated with the product 7T„: 



J? := 7 + mj. 



(1.24) 



1.5 Decompositions of SL(2,R) and the continuum regime 

Our next task is to define a useful continuum limit for the product 7T„. We 
remark that, for k = 1, the disordered quantum system introduced earlier 
expresses the matrix M as the product of a rotation matrix, a diagonal matrix 
and an upper triangular matrix: 

By applying the Gram-Schmidt orthonormalisation algorithm to the columns 
of M, it is easily shown that every element of SL(2,R) may be written in this 
form; this is a particular instance of the Iwasawa decomposition of a semi- 
simple Lie group into compact, Abelian and nilpotent subgroups [36]. Other 
familiar decompositions 

M = T 1 (t 1 )T 2 (t 2 )T 3 {t 3 ) (1.26) 

into one-parameter subgroups are the Gauss decomposition, where T\ is lower 
triangular, T 2 diagonal and T 3 upper triangular, and the Cartan decomposi- 
tion, where T\ and T% are rotation matrices, and T 2 is diagonal. Later in the 
paper, we shall work exclusively with the Iwasawa decomposition. Neverthe- 
less, there is some merit in explaining in some generality how a decomposition 
such as (1-26) determines a certain continuum limit; see for instance [37] for a 
study of the "Brownian bridge" limit of a product of random matrices where 
the Cartan decomposition is put to good use. 

Let us assume, then, that, for i = 1,2,3, the matrix Tj(tj) is a one- 
parameter subgroup such that 1^(0) is the identity. We choose as represen- 
tation space the set of functions on the projective space. Then the operator 
■%{ti) acting on functions via 

3r i {t i )f{z):=[foTr\t i )]{z) (1.27) 
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is a representation of the subgroup. Its infinitesimal generator is 



% := lim 

ti-s-0 



T = 5i(z) d? 



(1.28) 



where gi is a polynomial in z of degree at most 2 whose precise form depends 
on the details of the subgroup. 

In our context, the parameters U are random variables; we use the notation 

t7:=E(*i), Sti=ti-U, and Dij = E (5U 5tj) , (1.29) 

so that we have, in particular, 

ti=U + Sti. (1.30) 

It will be convenient to collect the means of the parameters into a vector and 
the covariances into a matrix: 



and 



D n D12 Dis\ 
D21 D22 D23 I 
, -D31 D 32 D 33 J 



(1.31) 



The covariance matrix is symmetric and non-negative; its square root er is 
therefore well-defined. 

The continuum regime is defined as the scaling regime where all the ex- 
pected values and the covariances are simultaneously small. The integral equa- 
tion (1.17) for the unknown density / may then be approximated by expand- 
ing the integral on the left-hand side in powers of the random parameters, 
neglecting terms of order greater than two, and taking expectations. The full 
derivation will be presented in Sec. 2.1. This yields the following first-order 
differential equation for /: 



_d 

dz 



/(*) 



- v{z) f(z) = j. 



(1.32) 



The local variance a 2 (z) and the local velocity v(z) appearing in this equation 
are given respectively by 



and 



a (z) := |erg(z) 



where c is the "correlation vector" given by 



c := 




(1.33) 
(1.34) 

(1.35) 



and g(z) is the vector- valued function with ith component g%{z). We shall refer 
to (1.32) as the integrated Fokker-Planck equation associated with the product 
of matrices. 
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1.6 A stochastic differential equation 

The continuum regime introduced above also makes sense in the non-stationary 
case. The continuum limit of (1.20) takes the form of the usual Fokker-Planck 
equation with z-dcpendent local velocity and local variance coefficients. An 
equivalent, more striking expression of this passage from the discrete to the 
continuous is the Stratonovich stochastic differential equation 

^ = fi ■ g(z(x)) - \ c • [g(z(x)) x g'(z(x))} + g(z(x)) ■ V (x), (1.36) 

where the continuous variable x replaces the discrete variable n, and the con- 
tinuous process z(x) replaces the discrete process {z n } defined by the random 
recursion (1.18). Here, the ith. component of the vector r\ is a white noise 
mix) associated with the component Ti(U) in the decomposition (1.26), the 
correlations between the components are given by 

E (7]i(x) r]j(x')) = D i:j 5(x - x') (1 < i, j < 3) (1.37) 

and, as before, the details of the decomposition (Iwasawa, Gauss, Cartan, etc.) 
are encoded in the vector g(z). 



1.7 The Iwasawa decomposition and the Hilbert transform 

Let us now restrict our attention to the particular case where the Iwasawa 
decomposition is used. For the ordering t\ — a, t 2 = w and t 3 = m, we have 




(1.38) 



We shall write D aa instead of Du, D aw instead of Di 2 , and so on. 

A "frontal attack" on the problem of computing the Lyapunov exponent 
would consist of finding the invariant density, and then evaluating the multiple 
integral in the Furstenberg formula (1.3). In the continuum regime, the right- 
hand side of this formula for the Lyapunov exponent reduces to a single integral 
with respect to the projective parameter z: 

/oo 
zf(z)dz 
-oo 

1 f°° d 

+ 2 Daa J ^[(i + ^/M] dz, (1.39) 

where the integrals on the right-hand side are Cauchy principal value inte- 
grals. Although this is a considerable simplification, the task of evaluating the 
characteristic exponent in terms of familiar functions is still quite daunting in 
the general case where the matrices M n do not commute among themselves. 
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In this work wc shall follow a different approach, in which a central part is 
played by the Hilbert transform F of the density 



F{V) ■= / 



m 

y- z 



dz, 



(1.40) 



where y is a complex variable in the lower half-plane. Either the above trans- 
form — whose relevance can be traced back to Dyson's seminal paper [5] — or 
its primitive 



J — c 



My - z ) f( z ) dz 



(1.41) 



have already been instrumental in the derivation of exact solutions of the 
Furstcnbcrg equation (1.5) in many situations [3,9-16]. 

We proceed to explain how the introduction of the Hilbert transform facil- 
itates our task in the continuum limit. Dividing the integrated Fokker-Planck 
equation (1.32) by (y— z) and integrating over z yields the following equation 
for F: 



d 

dy 



° 2 {y) 



F(y) 



-v(y)F(y) = n + ^(l + y 2 )+ay-w. 



(1.42) 



This equation for the Hilbert transform thus has the same homogeneous part 
as the integrated Fokker-Planck equation for the density; the crucial point is 
that the quantity of interest — the characteristic exponent — now appears on 
the right-hand side. 



1.8 Explicit form of the key equation 

To proceed, it will be useful to work with a fully explicit form of (1.42). For 
the Iwasawa decomposition, the local velocity and variance coefficients are 
polynomials in the Riccati variable z with respective degrees 3 and 4: 

v(z) = -a(z 2 + 1) + 2wz + u 
+ D aa z(z 2 + 1) + 2D ww z 
4-D tl l( ; z i 2 l D au z -j- *2D WU ^ 
a 2 (z) = D aa (z 2 + l) 2 + 4ZW 2 + D uu 

- 4D aw z(z 2 + 1) - 2£ Q „(z 2 + 1) + 4D wu z. (1.43) 

Hence the differential equation (1.42) for the Hilbert transform F(y) is equiv- 
alent to 

Q(y)F'(y) + R(y)F(y) = S(y) + 2f2, (1.44) 
where the coefficients 



Q(y) = a 2 (y), R(y) = 2 [a(y)a'(y) - v(y)} , 



(1.45) 
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and S(y) are polynomials in the complex variable y of degree 4, 3 and 2 
respectively. They are given explicitly by the formulae 

Q(y) = D aa (y 2 + l) 2 + 4D ww y 2 + D uu 

- 4D aw y{y 2 + 1) - 2D au (y 2 + 1) + 4D wu y (1.46) 

and 

R(y) = Ro(y) + Mv), S(v) = So{v) + S 2 (y), (1.47) 

where the polynomials with subscript are linear in the mean Iwasawa pa- 
rameters, and those with subscript 2 are linear in their covariances: 

Ro(y) = 2(ay 2 - 2wy + a-u), 
S (y) = 2(ay-w), 

R 2 (y) = 2(D aa y 3 - 2D aw y 2 + (D aa + 2D ww )y - 2D aw ), 

S 2 (y) - D aa (y 2 + 1). (1.48) 

1.9 Outline of the solution method. Halperin's example 

The differential equation (1.44) can be solved in two steps. First, the inte- 
grating factor — the solution of the corresponding homogeneous equation — is 
given by 

H(y) ^ cxp {-Iw) dy )- (L49) 

Second, the solution of the full, inhomogeneous equation can be obtained by 
"varying the constant" : by setting 

F(y) = K(y)H(y), (1.50) 

Equation (1.44) is reduced to 

_ S{y) + 2f2 

{v) ~ Q(y)H(y) ■ (L51) 

We thus obtain 

The function F{y) that we seek has two distinguishing properties: firstly, 
it is analytic in the lower half-plane. Secondly, since it is the transform of a 
probability density, it has the asymptotic behaviour 

F(y)nl/y as |y| ^ oo (1.53) 

along any ray contained in the lower half-plane. 1 These two properties are 
sufficient to determine both the characteristic exponent £2 and the interval of 
integration in (1.52), i.e., the full function F{y). 



1 In this work the notation A(x) B(x) as x — ¥ xq means lim x _» xo (A(x)/B(x)) = 1. 
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Let us illustrate this all-important point with a concrete example. Consider 
a product of matrices of the form (1.25) where only the u n are random. This 
product corresponds to a disordered quantum system with a potential made 
up of delta scatterers whose strengths are random, and the continuum regime 
is effectively the limit studied by Halperin [33] . The only non-zero parameters 
are a, W, u and D uu . Hence 

Q(y) = D uu (1.54) 
is a constant. Equation (1.49) yields 

H(y) = c p (y)/ D ^ j (1.55) 

where P(y) is the polynomial 

P{V) = -^ay 3 + 2wy 2 + 2(u-a)y, (1.56) 

whilst Equation (1.51) reduces to 

K'{y) = -^-(ay -w + n)e~ p ^! D ^ . (1.57) 

Assume for dcfinitcness that a > 0. Then H exhibits exponential growth 
at infinity along any ray contained in the unshaded sectors of the complex 
plane depicted in Figure 1.1. Since F is analytic in the lower half-plane, Equa- 
tion (1.50) and the asymptotic condition (1.53) together imply that K must 
decay to zero along any ray in the lower unshaded sectors. Take for instance 
the piecewise linear path shown in Figure 1.1: it originates at — oo — iO and 
follows the real axis in one of the unshaded sectors until it reaches the origin, 
where it is refracted and continues to infinity along the ray argy = — n/3 in 
an another unshaded sector. K vanishes at the extremities of this integration 
path; using (1.57) and the Fundamental Theorem of Calculus, we deduce 



W3 r 



/e ' oo 
(ay - w + Q) c - p( *y)/ D ^ dy = 0. 
-OO 

The linear change of variable 



(1.58) 



Pt + w 

y = _ , (1-59) 



a 



where 



V3 and 



IS 

(see (3.8)), brings the contour integral into the following form: 



au-a z )=^ (1.60) 



L 



{pt+ f2)e t3 / 3 ~ xt dt = 0. (1.61) 
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Imy 




Rcy 



Fig. 1.1 Path of integration for monolithic distance disorder. The integrating factor H 
exhibits exponential growth at infinity in the unshaded sectors, and exponential decay in 
the shaded ones. 



The integrals involved can be expressed in terms of the Airy function [38, 
(9.5.4), p. 196] 



t 3 /3-xt 



and of its derivative with respect to x. The end result is 

Q = pG{x\ 
where the scaling function G(x) is 



G{x) = c 



- n -2V3 



Ai'(e- 2i7r / 3 x) _ Ai'(x) + iBi'(x) 
Ai( c -^/3 x ) ~ Ai(ic) + iBi(ic) ' 



(1.62) 
(1.63) 

(1.64) 



This formula will be analysed in Sec. 4.1. 

In the general case, the complex zeros of the polynomial Q(y) — which is 
nothing but the complexified diffusion coefficient — will be at the heart of the 
subsequent analysis. It is indeed clear from (1.49) that the auxiliary function 
H(y) is singular at the zeros of the polynomial Q(y), and that the multiplici- 
ties of the zeros dictate the nature of the singularities. For instance, a simple 
zero at y\ yields the power-law singularity H(y) <~ (y — y\) a , a double zero 
at 2/1 yields the essential singularity H(y) <~ (y — yi) a cxp(6/(y — yi)), and 
so on. 2 It will be advantageous to take the view that the polynomial Q(y) 



2 In this work the notation A(x) ~ B(x) as x — !• xo is weaker than A{x) ~ B(x). It means 
that A(x) and B(x) become proportional as x — > xo, up to an unspecified factor, which may 
either be constant or vary much more slowly than A(x) or B(x). 
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always has four zeros, counted with their multiplicities — possibly including 
a zero at infinity. Thus, for Halperin's example, we would say that Q has a 
quadruple zero at infinity. This convention is justified by the fact that the 
Riccati variable z is a projective coordinate. Table 1.1 summarizes the rela- 
tionships we have found between these zeros, the types of disorder, and the 
special functions that describe the scaling form of the characteristic exponent 
in the continuum regime. 



Pattern of zeros 


Type of disorder 


Special function 


Section 


1 quadruple 


scalar 


Airy 


4.1 


2 double 


supcrsymmctric 


Bessel (real index) 


4.2 


2 double 


distance 


Bcsscl (imaginary index) 


4.3 


1 double + 2 simple 


general potential 


Whittakcr 


5 


4 simple 


independent (zero mean) 


elliptic 


6 


4 simple 


most general 


hypcrgcomctric 


7 



Table 1.1 List of the examples worked out in this paper, giving the correspondence between 
the pattern of zeros of the polynomial Q(y), the type of disorder, and the type of special 
function entering the scaling form of the characteristic exponent Q. The last column indicates 
the section where the details arc carried out. 



2 Basic tools and derivation of the key formulae 

In this section we present a detailed derivation of the key formulae upon which 
our study relies, namely of Equations (1.32), (1.39) and (1-42). 

2.1 The equation for the density in the continuum regime 
We use the notation of Sec. 1.5 and set 

[M-\z) 

Z(h,t2,t 3 ) := dt.f{t) = P{M- 1 {z))-P(z), (2.1) 

J z 

where P is a primitive of /. Our aim is to determine the quadratic Taylor 
polynomial of £ when the origin is the point of expansion. An elegant way of 
doing this is to introduce the representation 2T of SL(2,R) which maps the 
matrix 

M = T 1 (t 1 )T 2 (t 2 )T 3 (t 3 ) (2.2) 
to the operator £?m defined by 

&m:= &l{H)o £T 2 {t 2 )o *r 3 {t 3 ). (2.3) 

Then 

£(*i.*2,t3) - (z) - {STiP) (z) (2.4) 
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where I is the identity matrix. The operator 3?m may be expressed in terms 
of the subgroup generators as 



2?m = e* 1 ^ 1 e* 2 ® 2 e' 3 ® 3 . 



(2.5) 



Reporting this in (2.4) and expanding the exponential, we obtain 

£(*l,*2,t3) 



+tl*2 ^1 ^2 + £l^3 ^1 ^3 + *2*3 ^2 ^3 + ' ' ■ 



P(z). (2.6) 



The identity 3>iP{z) — f(z)@iZ then leads to the following formulae for the 
partial derivatives of £: 



dU 



tl=t 2 =t3=0 



and 



9 2 C 



dtidtj 



%\f{z)^z\ 



if « < 3, 



ti=t 2 =t 3 =0 



\f{z)9jz] + f(z)[9i,9j]z if i> j, 



(2.7) 



(2.8) 



; -] is the commutator of ^ and Replacing £ in (1.17) by its 



where 

quadratic Taylor polynomial leads to (1.32). 



2.2 Calculation of the invariant measure 

Although our approach does not require explicit knowledge of the invariant 
density /, its calculation is often of independent interest, and so we discuss it 
here. As a by-product, we shall obtain information on the large- z behaviour 
of f(z) that will prove useful in the derivation of the equation for its Hilbert 
transform. 

The solutions of the Fokker-Planck equation form a two-parameter family: 
the current j is one of these parameters; the other is a constant of integration. 
We now explain how these two parameters must be chosen in order to yield 
the invariant measure. It will be helpful to work with the angle variable 

ip = 2 arctan z € [— x, x] . (2.9) 

Let us indicate very briefly the equations appropriate for this choice. We shall 
use the subscript a to distinguish the density f a of the angular variable ip 
from the density / of the projective variable z, and similarly for the angular 
counterparts of other functions such as a 2 , v, g, etc. We can express the 
relationship between the densities as 



(l + z 2 )f(z) = 2f a (p). 



(2.10) 
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In the continuum regime, it is easy to see by substitution in (1.32) that the 
density f a satisfies 



d 

dip 



U<p) 



v a (<p)f a (ip) 



(2.11) 



where the local variance a\ and the local velocity v a are given by formulae 
analogous to (1.33)-(1.34): 



(2.12) 



and 



««(¥>) ~ g {s'aif) ■ <r 2 ga&) - c • [g a {ip) X g^)] } - (1 ■ g a {ip), (2.13) 
where 

8a(<p) = -2sinip . (2.14) 

y — i — cos ip J 

In the continuum regime, the angle variable is therefore a diffusion process on 
the interval [— it, it] with infinitesimal generator 



d 2 , . . d 



2 d<^ 2 



acting on the space of twice-differentiable 27r-periodic functions. 
Set 

JM:=2r^|d, 

— 7r ff aW 

Equation (2.11) may be integrated to yield 

2e J (v) 



/.fa) 



C + j 



J « dt 



(2.15) 



(2.16) 



(2.17) 



The two parameters C and j in this expression are determined by requiring 
that f a be normalised and 27r-periodic. The periodicity condition yields 



1 - e J « 



J — IT 



d<p>. 



(2.18) 



There are two cases to consider. 
• The first case arises when 



J— 7T "o 



(*) 



di = 0. 



(2.19) 



Then the periodicity condition (2.18) implies j = 0. The general situation 
where the stationary probability current vanishes will be investigated in 
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Sec. 7.3, while we shall encounter an interesting special case in Sec. 8. For 
the time being, let us observe that the condition (2.19) is equivalent to 



o»(z) 



0. 



The density is then given by 



fa{<f) = 2C 



e J(<p) 



(2.20) 



(2.21) 



and C is the normalisation constant. 
• The second case corresponds to J(ir) ^ 0. Then the periodicity condition 
(2.18) expresses C in terms of j, the density is given by 



fa(<p) = 



2.7 



[IX pip 

e J W / e~ J ^dt + / , 

J if J — 7T 



1-e'M alb) 
and j plays the role of a normalisation constant. 



(2.22) 



2.3 The Rice formula 

By setting ip — ±7r in (2.11), we readily obtain 

iD aa f a {±7r) + 2(a + 2D aw )f a {±7r)=j. (2.23) 

When D aa > 0, the periodicity of f a therefore implies the periodicity of (all) 
its derivative(s). By expressing this periodicity in terms of the density / of the 
Riccati variable z, we find 

lim z 2 f(z) = lim z 2 f(z) (2.24) 

z— > — OO z— >OQ 



and 

~ "f [(l + z 2 )f(z)] = limz 2 -^ 

z^-<x> dz 1 J z^oo dz 



Furthermore, (2.23) becomes 

D aa lim z 2 ^-[(l + z 2 )f(z)]+(a + 2D aw ) lim z 2 f(z)=j. (2.26) 

|z| — ^oo QZ |z| — >-oo 

In particular, the case where a is non-random (i.e., D aa = 0) yields 

lim z 2 f(z) = =. (2.27) 
|z|->oo a 

We thus recover a well-known relationship between the integrated density of 
states and the decay of the stationary density at infinity [32] . The above iden- 
tity can be viewed as a special case of a formula due to Rice for the density 
of zeros of a continuous process [39]. We shall henceforth refer to (2.26) as the 
Rice formula. 
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2.4 The Lyapunov exponent in the continuum regime 
Following [21, Sec. 3.4], we write 

so that the Furstenberg formula (1.3) for the Lyapunov exponent may be 
expressed in the form (see the bottom of p. 441 in [21]) 

i / r°° d r B ' 1(z) \ 

7 = -™+-E(/ dzln(l + z 2 )-^ d*/(t)J. (2.29) 
By using integration by parts in the integral over z, we obtain 

/ rco rB-\z) \ 

7 = -1J - E (/ d* Y — 5 jf dtf(t) j . (2.30) 

This formula is valid whether or not one considers the continuum regime; we 
shall now show how it simplifies in that regime. 
Set 

x ~L*t+?L dtm - (23i) 

Just as, in Sec. 2.1, we expanded the integral £ in powers of a, w and u, we 
can expand the integral \ in powers of w and u. After taking expectations, we 
find 

/CO ^ 
dz 2 h(w,u,D ww ,D wu ,D uu ;z) , (2.32) 
-oo 1 + Z 

where 

h := {w9 w z + u£> u z + ^D wu \® w , S> u ] z^j f(z) 

+ ~{d ww S> w [f{z)3> w z] + D WU S> W [f(z)3> u z] 

+ D wu $> u [f(z)S> w z] + D UU 9 U [f{z)® u z]}. (2.33) 

We then use the integrated Fokker-Planck equation (1.32) to express h in 
the equivalent form 

h = j -af(z)@ a z-^@ a {f(z) {D aa @ a z + 2D aw @ w z + 2£> au # ti z]}. (2.34) 

Here, the important observation is the occurrence of the operator & a in all 
but the first term on the right-hand side. Since 

® a = (l + z 2 )-^-, (2.35) 
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the awkward denominator in the integrand of (2.30) will be removed. Report- 
ing this in (2.32), we find 



7 = — w + a 



+ 



f 

J — C 



zf(z) dz 



(z) + D 9u(z) 



dz,(2.36) 



where the gi(z) = %z, i £ {a, w, u}, are the components of the vector (1.38). 

This formula may be simplified somewhat by performing an integration by 
parts involving the last two terms in the square brackets; this yields (1.39). The 
derivation of (1-42) for the Hilbert transform is now reasonably straightforward 
but somewhat lengthy (see Appendix A). 



3 The weak-disorder regime 

We describe in this section the systematic expansion of the characteristic ex- 
ponent in powers of the covariances. Although the expansion is divergent, the 
knowledge of the first few terms has its own interest. Besides this, it will also 
provide a useful means of checking the correctness of the scaling forms that 
we shall determine later on. 



3.1 No disorder 

We start by considering the case where there is no disorder. In this situation, 
all the matrices M n are equal to the constant matrix Mq corresponding to a, 
w and u. The latter matrix is of the form (1.7) with 

mn = cos a e w , m\i — cosa e w u — sina e~ w , 

TO21 = sinae™ , m,22 = cosSc - " 1 + sinae^TL (3.1) 

Let us introduce the variable /x such that 

tr Mo = sin a e w u + 2 cos a cosh w = 2 cosh /i. (3.2) 

The eigenvalues of M are therefore e ±M . The recursion (1.18) reduces to the 
deterministic Moebius transformation 

z n +i=Mo(z n ) (3.3) 

whose fixed points are 

cosa-e-C^) 

z± = — . 3.4) 

sin a 

Let us assume for definiteness that \x has a non-zero real part, and choose fi 
so that 

Reit>0. (3.5) 
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If |trMo| < 2, so that (3.2) yields imaginary values of /x, we supplement /j, 
with an infinitesimal positive real part, so that (3.5) still holds. This condition 
ensures that z + (resp. z_) is the stable (resp. unstable) fixed point of the 
mapping Mo- We thus obtain the simple result 

/„(*) = S(z - z+) (3.6) 

for the invariant density of the Riccati variable, and the expected expression 

= H (3.7) 

for the characteristic exponent. 

In the continuum regime where a, w and u are small, the above expres- 
sions (3.2) and (3.4) respectively simplify to 



w 2 + au — a 2 



(3.6 



and _ _ _ 

w ± u a — u . . 

z± = — =^ = = . 3.9 

a to =p /J, 

Throughout the following, \x will be a notation for the right-hand side of (3.8). 

Let us now examine the general formalism in the present non-random sit- 
uation. The only non-zero polynomials are Ro{y) and So(y)- Equation (1.44) 
therefore loses its differential character. Using the subscript to indicate that 
there is no disorder, we obtain 

The polynomial Ro(y) factorizes as 

Ro(y) = 2a(y-z+)(y-z-), (3.11) 

where z± are the fixed points (3.9) of the mapping (3.3) in the continuum limit. 
The natural condition that F (y) have no pole at the unstable zero (y = zJ) 
allows one to recover the results (3.6) and (3.7). 



3.2 Systematic weak-disorder expansion 

The above line of thought can be pursued and yields a systematic weak- 
disorder expansion to all orders in the continuum regime. 

Let us look for a perturbative solution to (1.44) of the form 

F(y) = F {y) + F 2 {y) + F 4 (y) + ■ ■ ■ , [2 = Q u + f2 2 + Q A + ■ ■ ■ (3.12) 

The functions i r 2fc(y) and the corresponding contributions Q 2 k to the weak- 
disorder expansion of the complex characteristic exponent can be obtained 
recursively. 
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3.2.1 Second order 

The equation for F 2 (y) reads 

Ro(y)F 2 (y) = -Q(y)F^(y) - R 2 (y)F (y) + S 2 (y) + 2J7 2 . (3.13) 
The condition that F 2 (y) have no pole at the unstable zero (y — z_) yields 

^2 = l -[Q{y-)F^y_) + R 2 (z_)F (z_) - S 2 (z-)] . (3.14) 

We thus obtain the following explicit expression for the second-order contri- 
bution Q 2 : 

-8n 2 f2 2 = {Aw 2 + u 2 ) D aa 
+ Aa(a - u) D ww 
+ a 2 D uu 

+ 4(— ufi — 2aw + wu) D aw 
+ 2(2wJ/Lt — 2w 2 — au) D au 

+ Aa(w - n) D wu . (3.15) 

We also obtain an explicit rational expression for F 2 (y), with a triple pole at 
the stable zero (y = z + ). 

It is worth pointing out that the result (3.15) coincides with the continuum 
limit of the following general expression for the second-order weak-disorder 
expansion of the complex characteristic exponent for arbitrary values of the 
mean variables a, w and u: 



—8 sin a sinh 2 /x Q 2 = sina(u 2 e 2w + 4sinh 2 w) D aa 
+ 4 sin 2 a(sin a — u cos a) D ww 
+ sm 3 ae 2w D uu 

— 8(cosh(w — /x) — cosa)(l — cosac~ M '~ Ai ) D aw 

— 4sina7c"'~ AI (cosh(w — n) — cos a) D au 

— 4 sin 2 a(cos a - e w ~^) D wu . (3.16) 

This result can be derived by standard methods [3,24]. 

3.2.2 Fourth order 

The equation for F±{y) reads 

R (y)F 4 (y) = -Q(y)F! 2 {y) - R 2 (y)F 2 (y) + 2f2 4 . (3.17) 
The condition that F 4 (y) have no pole at the unstable zero (y = z_) yields 

^4 = \ [Q(z-)F 2 (z-) + R 2 (z_)F 2 (z-)] . (3.18) 
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We thus obtain the lengthy expression 

-128m 5 ^4 = (4uJ 2 + u 2 )(165 2 - 165w + 4w 2 + 5u 2 )D 2 al 
+ 165(5 — u)(a 2 — au + Aw 2 )D 2 vw 
+ ^D 2 UU 

+ 16 [Auw{2a - u)n + 20a 2 w 2 

+a 2 u 2 - 205 w 2 u - au 3 + Aw 2 u 2 ] D 2 
+ 45 [ - 8w(2a - u)fi + m 2 u + 20a w 2 - 3au 2 - 8uw 2 ] D 2 al 
+ 165 2 (-4w/i + a 2 -au + Aw 2 )D 2 wu 
+ 8(205 2 w 2 - a 2 u 2 + au 3 - 20auw 2 + 6u 2 w 2 )D aa D ww 
+ 2(8a 3 u + I2a 2 w 2 - 3a 2 u 2 + 8w 4 )D aa D uu 
+ 8(25 - u) [ - 2u(25 - u)fi 



) 2 

aw 



-w(8a 2 - 8au + 5u 2 + 12w 2 )] D aa D aw 



+ 4 [4uJ(25 - u) 2 M - 8o^u - 325 2 w 2 + Aa 2 u 2 - au 3 

+20auw 2 - 8TZJ 4 - 6u 2 w 2 ]D aa D au 
+ 8[- 2(25- u)(2w 2 + au)(j, 

+w(8a 3 + 12aw 2 - 3au 2 - 8uw 2 )} D aa D wu 
+ 85 2 (-5 2 + au + 6w 2 )D ww D uu 

+ 32w{—2wuu — 65 3 + 9a 2 u — Aaw 2 — 3au 2 + 2uw 2 )D ww D aw 

+ 165 [4(5 — u)wu + a 2 u — 10a w 2 — au 2 + Au w 2 ~\ D ww D au 

+ 32aw(—2w(i + 3a 2 — 3au + 2w 2 )D ww D wu 

+ 8a [2(-25 2 + au + 2w 2 )u - w{2a 2 + 3au + 8w 2 )] D uu D aw 

+ 45 2 (4W/x — 4a 2 -51- &w 2 )D uu D au 

+ 8a 3 (hw-2u)D uu D wu 

+ 16 [2(25- u)(au- 2w 2 )u 

+w(8a — Aau + au + 2aw — Auw )\D aw D au 
+ 32{AW 3 f i - 45* + ha 3 u - 25 2 uJ 2 - 5 2 w 2 - Aw^D^D^ 
+ 165[2(25 2 - au + 2w 2 )u 

+w(-6a 2 + au- Aw 2 )] D au D wu . (3.19) 

We also get an explicit rational expression for F\(y), with a pole of order 5 at 
the stable zero (y = z+). 

3.2.3 Generic structure 

The structure of the weak-disorder expansion clearly appears from the first 
two orders studied above. At a generic (even) order 2k for k > 2, the equation 
for F 2 k(y) reads 

Ro(y)F 2 k(y) = -Q{y)F! ik _ 2 {y) - R 2 {y)F 2k _ 2 (y) + 2f2 2k . (3.20) 
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The condition that F 2 k{y) have no pole at the unstable zero (y = zJ) yields 
V 2k = \ [Q(z-)F^ k _ 2 (z_) + R 2 (z_)F 2k _ 2 (z_)] . (3.21) 

The structure of the expression thus obtained is clear from the recursive nature 
of the problem. The denominator of f2 2k is an integer multiple of /x 3k ~ 1 . The 
numerator is a homogeneous polynomial in the six covariances of degree k. 
The explicit results obtained above at orders 2 and 4 suggest that the 

A ( fc + 5 ) ! (* 99) 

dk ~ Hd5T (3 - 22) 

different possible monomials are all present in the end result. We have d\ = 6, 
d 2 = 21, g?3 = 56, c?4 = 126, and so on. The coefficient of each monomial 
is a homogeneous polynomial in a, w, u and fi, of degree 2k with integer 
coefficients. The degree in /i can be reduced to one by means of (3.8). Finally, 
Fik{y) is a rational function, whose only pole is a multiple pole of order 2k + 1 
at the stable zero (y = z + ). 



4 Monolithic disorder 

In this section we deal with the situations where only one of the three variables 
a n , w n or u n is random. We refer to these three special cases as monolithic 
disorder. 



4.1 Scalar disorder (only the u n are random) 

We began our study of this case in the introduction, where we derived the 
following exact formula for the characteristic exponent (see (1.63) and (1.64)): 

Q = pG(x), (4.1) 

where the scaling function G(x) is 

_ W3 Ai'(c-2W3 x) _ Ai'QO+iBi'Cs) 
{ ) M(e-**/ 3 x) Ai(x)+iBi(x) ' 1 ' 



with 



— 2 n \ 1/3 2 

a z D uu \ n 2 

and x = (4.3) 



The argument x of the Airy functions involves /U and (3. The first of these 
parameters is a measure of the distance to the band edge (we have /i 2 < in 
the band and fi 2 > in the gap), while the second one demonstrates that the 
effective disorder strength is the product a 2 D uu . The precise correspondence 
with the quantum-mechanical problem is discussed in Sec. 9.3. 
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The Wronskian identity [38, (9.2.7), p. 194], [40, (10.4.10), p. 446] 



yields in particular 



Ai'(x)Bi(x) - Ai(a;)Bi'(x) = - (4.4) 

7T 



lmG(x) = , ... ,„ — —. .„. . (4.5) 

V ' 7T [Al(x) 2 + Bl(x) 2 ] V ' 



The above Airy scaling is characteristic of a generic band edge singularity. 
It has been met in a variety of situations, including the problem of a white 
noise scalar potential [32,33,10,41] and the weak-disorder regime of the tight- 
binding Anderson model with diagonal disorder near its band edges [42] . 

There remains to extract from these formulae the behaviour of the charac- 
teristic exponent in various interesting limits. The differential equation obeyed 
by the Airy functions [38, (9.2.1), p. 194], [40, (10.4.1), p. 446] 

Ai"(x)-xAi(x) =0 (4.6) 

translates into the following Riccati equation 

G(x) 2 + G'(x) = x (4.7) 

for the scaling function G(x). Solving the above equation iteratively yields the 
asymptotic expansion 

15 15 

G(x) = x 1 ! 2 - - r-r^r - —— H as x ->• oo, (4.8) 

v 1 Ax 32x 5 / 2 64x 4 ' v ; 

which corresponds to the weak-disorder expansion 

_ a*D uu ^Dj u 

8/x 2 128^ 512^ + as ^»^ u - ^> 

The first two non-trivial terms are in agreement with (3.15) and (3.19). 

It should be remarked, however, that the asymptotic expansion (4.8) is 
divergent. In particular, it fails to capture the exponentially small imaginary 
part of the form 

ImG(x) ~ exp (-^ 3/2 ) (4-10) 

of the scaling function for large positive values of x, which describes the tail 
of the density of states deep in the gap. 

We are led to draw from this explicit example the conclusion that the 
weak-disorder expansion is generically only a divergent asymptotic expansion. 

Right at the band edge (x — 0), we have 

_ 2W3 Ai^o) _ W3 r(2/3) ( m 2 D uu \ 1/3 

Ai(o)^™ r(i/3)\ 2 ) ' (4 - iij 



25 



and hence 

r(2/3) f 3a 2 D uu \ 1/3 . V3r(2/3) f 3a 2 D uu \ 1/3 

^=2Tm{-^) and j = 2^T{im{-^) • (4 - 12) 

To close, let us consider the case where a = w = 0. In this situation, the 
matrices 

M n = T 3 (u n )= (I (4-13) 

belong to the nilpotent subgroup of SL(2, R), and hence commute among them- 
selves. We thus have i7„ = T 3 (U n ) with U n = u\ H + u n « nl, and so 

17 = 0. (4.14) 

This result is recovered from (4.1), since /? = 0. 



4.2 Supersymmetric disorder (only the mj„ are random) 

We continue our study with the case where only the supersymmetric vari- 
ables w n are random. The non-zero parameters are therefore a, w, u and 
D ww . The polynomial 

Q(y) = 4D ww y 2 (4.15) 

has a double zero at the origin and a double zero at infinity. Equation (1.49) 
yields 

H(y)=y v - 1 e*W, (4.16) 



with 

w 



(4.17) 



and 

$(iA = - 

2D 



*(y) = -7^[ay + ^-^\. (4.18) 



Equation (1.51) reads 

2D ww K'(y) = (ay-w + Q) y-^ 1 e"*^. (4.19) 

This function K{y) must vanish both as \y\ ~ > and \y\ — >• oo in the directions 
where <P(y) diverges. Assuming u > a > for definiteness, and changing the 
sign of the integration variable y for convenience, we obtain 

/>oo 

/ (ay + w-f2)y- v - 1 c i ' ( -y'>dy = Q. (4.20) 
Jo 

By setting t — 1/y, the integrals involved in the latter condition can be reduced 
to integrals of the form [43, vol. I, (29), p. 146] 



Jo 



//,W 2 

t „-i e - at -b/t dt = 2 I _\ Kv (2Vab), (4.21) 
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where K v is the modified Bessel function. Using the identities [44, vol. II, 
(25), (26), p. 79] 



xK v _\{x) + vK v {x) = xK v+ i(x) - vK v {x) = -xK'„(x), 
we eventually find 



ft = D ww G(x) 



where 



and (see (4.17)) 



G{x) 



K v {x) 



w 



D„ 



x = 



yj a(u — a) 



D„ 



(4.22) 
(4.23) 

(4.24) 
(4.25) 



The characteristic exponent Q is real and the stationary current j vanishes as 
long as the argument x is real, i.e., a(u — a) > 0. 

In the opposite situation, i.e., a(u— a) < 0, it is more convenient to consider 
the real variable 



c 



y/ 'a(a — u) 



D„ 



The result (4.23) becomes 



with 



Gi(C) = -C 



Q — D ww Gi(C), 



J v {Q+iN v {Q' 



(4.26) 



(4.27) 



(4.28) 



where J„ and N v are Bessel functions. The Wronskian identity [44, vol. II, 
(28), p. 79] 

^(0^(0-^(0^(0 = 4 ( 4 - 29 ) 



yields in particular 



ImGi(C) = 



tt(JAC) 2 + N v (0 2 Y 



(4.30) 



The above Bessel scaling functions have already been met in several cir- 
cumstances. The forms (4.28) and (4.30) enter the analysis of classical diffusion 
in a one-dimensional random force field [45,46], whereas the form (4.24) shows 
up in disordered supersymmetric quantum mechanics (see Sec. 9.4) and in the 
quantum Ising chain in a disordered transverse magnetic field [17]. 

The differential equation satisfied by the modified Bessel function [44, 
vol. II, (11), p. 5] 



x 2 K'l{x) + xK' v (x) - (x 2 + v 2 )K v {x) = 



(4.31) 
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translates into the following Riccati equation 

G(x) 2 - xG'(x) = x 2 + v 2 = (4.32) 

WW 

for the scaling function G(x). We use this equation to investigate the weak- 
disorder regime where 

A(x) = yjx 2 + v 2 = (4.33) 
is large. Solving (4.32) iteratively yields the asymptotic expansion 



This corresponds to the weak-disorder expansion 

a(u-a)D ww a(u-a)(a 2 -au + Aw 2 )D 2 ww 

v + +••• (4.35) 

as D ww — > 0, in agreement with (3.15) and (3.19). 

In the opposite strong-disorder regime, the expansion at small x 

C(*) = M + 2(TTH) + --- (436) 

translates into 

^Uww 

Let us turn to the case w = 0, where the index v vanishes. This situation, 
which will also be met for distance disorder with u = 2a, corresponds to a 
critical point. In the present case we have 

x = (4.38) 



D 



W W 



while the scaling law (4.23) becomes 



« = (4-39) 



K (x) 

The weak-disorder expansion (4.35) simplifies to 

n = fl+ Dww_D^ + ^ (440) 

In the opposite regime of a strong disorder, corresponding to x — > 0, the loga- 
rithmic singularity Kq(x) ~ ln(2/x) — C translates into the singular behavior 
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where C denotes Euler's constant. 

To close, let us consider the case where a = u = 0. In this situation, the 
matrices 



M n =T 2 (w n ) = 




(4.42) 



belong to the Abelian (diagonal) subgroup of SL(2,M), and hence commute 

among themselves. We thus have 77„ = T2(W n ) with W n = w\-\ \-w n w riw, 

and so 

Q = \w\. (4.43) 
This result is recovered from (4.36), since x = 0. 



4.3 Distance disorder (only the a n are random) 

We close this section with the case where only the distance variables a n are 
random. The non-zero parameters are therefore a, w, u and D aa . The poly- 
nomial 

Q(y) = D aa (y 2 + l) 2 (4.44) 
has two double zeros at y = ±i. Equation (1.49) yields 



1 - V 

l + iyj y 2 + l 



H(y)=(TT^.) Z^TT- (4-45) 



with _ _ 

Equation (1.51) reads 

D aQ ^'(l,) = (5(|,) + 2/2)^^ J ^ TT , (4.47) 

where S(y) is given by (1.48). 

The analysis is facilitated by working with the variable 

t = i±i» (4.48) 
1 - vy 

which maps the lower half-plane (Imy < 0) onto the complement of the unit 
circle (\t\ > 1), and the double zeros (y — i and y = — i) to zero and infinity 
respectively. The transformed function K(t) := K(y) satisfies the differential 
equation 



\D aa K'{t) = 



{Q — w — ia)t + {Q -w + ia)j + 2{Q - w + D aa ) 
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with <P(t) := <P(y), i.e., 
$(t) = 



4D n 



(2w + \u)t + (2w -iu)- + Aw 



(4.50) 



Let us assume w < for definiteness. The function K (t) must vanish as 
Ret — > — oo. The characteristic exponent fl is therefore determined by the 
condition 

K'{t)dt = 0, (4.51) 



where the integration contour C circles around the branch cut of the integrand 
along the negative real axis. 

An integration by parts, corresponding to the choice 



g(y) = 



y 



y 



l 



(4.52) 



of the gauge function (this trick will be explained in detail in Sec. 7; see 
Equation (7.11)), simplifies the above equation to 



(2w + m)t + (2w - iu) - - 4J? 



tiA-l e -*(t) dt = Q_ 



(4.53) 



The integrals involved in the latter expression are of the form 

J t v - 1 c at+b ' t &t = 2ni(^j ' I v (2^fab), (4.54) 

where I u is the modified Bessel function. This result, which bears a close 
resemblance with (4.21), can be easily derived from [47, (8.412.2), p. 954]. 
Using the identities [44, vol. II, (23), (24), p. 79] 



xl v+ \(x) + vl v {x) = xl„-i(x) - vl v {x) = xl' v (x), 
we eventually obtain 



f2 — D aa G(x), 



where 



and (see (4.46)) 



G(x) = x 



A 



u — 2a 



\jAw 2 + u 2 



(4.55) 
(4.56) 

(4.57) 
(4.58) 



ID 2D 

The present case is in some sense dual to the previous one. Both scaling func- 
tions (4.24) and (4.57) involve a modified Bessel function. In the present situ- 
ation of distance disorder the argument x is real and the index v = iA is imag- 
inary, whereas in the previous case of supersymmetric disorder (see (4.25)), x 
can be either real or imaginary and the index v is real. 
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The differential equation (4.31) is also satisfied by I v {x). We can study 
the weak-disorder regime by using the corresponding Riccati equation for the 
scaling function G{x): 

G{x) 2 + xG'{x) = x 2 — \ 2 = YpT- ( 4 - 59 ) 

act 

Set 

A{x) = Vx 2 - A 2 = (4.60) 
By solving (4.59) iteratively, we obtain the asymptotic expansion 

x 2 x 2 (x 2 +4:X 2 ) 

This corresponds to the weak-disorder expansion 
(4W 2 + u 2 )D aa 



(4wT 2 + u 2 )(4w 2 + Su 2 - 16au + 16a 2 )D 2 aa 



as D aa — > 0, in agreement with (3.15) and (3.19). 

In the strong-disorder regime, the small- a; expansion 

x 2 



G{x) = -iX +w -r r + ... (4.63) 



translates into 



.2a- u Aw 2 + u 2 /„ „^ 

J? = i — - — + — — + • • • as D aa -> oo. (4.64) 

I oL> aa 

Let us turn to the case where u = 2a, so that the index A vanishes. In this 
case we have 

x=-^-, (4.65) 
while the scaling law (4.56) becomes 

The weak-disorder expansion (4.62) simplifies to 

« = ^~ ~ +- as D aa -+ 0. (4.67) 

In the opposite regime of a strong disorder, corresponding to x — > 0, we obtain 
the estimate 

u 2 

Q w -p— as Daa -> oo. (4.68) 
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To close, let us consider the case where w = u = 0. In this situation, the 
matrices 

M n — Ti(a n ) — ( C ° S ~ Sm an ) (4.69) 
y sin a n cos a n J 

belong to the compact (rotation) subgroup of SL(2,R), and hence commute 

among themselves. We thus have 7T„ = T\{A n ) with A n = a\ H \-a n rts ria, 

and so 7 = 0, while j = ±a/n. The sign of j depends on how the limiting case 
is reached. Our prescription W — > 0~ yields 

Q = iTrj = ia. (4.70) 

This result is recovered from (4.63), since x = 0. 



5 General potential disorder 

In this section we deal with the case of a general potential disorder, whose 
supersymmetric and scalar parts 8w n and 5u n are correlated. The non-zero 
parameters are therefore a, w, u and the three covariances D WWl D wu and D u 
For further convenience we introduce the reduced quantities 



c =-====, ^Vl-c 2 , yo = -J——. ( 5 -!) 
The polynomial Q(y) then reads 

Q(y) =4 J D ww ( 2 / 2 + 2cy y + yo)> (5-2) 
where c is nothing but the correlation coefficient between 6w and Su. 



5.1 The fully correlated cases (c 2 = 1) 

Let us consider the fully correlated cases where c = ±1. The polynomial 

Q(y) = 4D ww (y + cy ) 2 (5.3) 

has a double zero on the real axis at y = —cy and a double zero at infinity. 
Equation (1.49) yields 

H(y) = (y + cy o r- 1 c^y\ (5.4) 

where 

$(y) = -— j— (ay+ - V (5.5) 
W 2D WW \ y y + cyoj V ' 



and 



f= — ^ , P = u - a + 2(D WW - w)cy - ay a . (5.6) 
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Equation (1.51) reads 

2D ww K\y) = (ay-W + Q){y + cy^- 1 e~*M. (5.7) 

Assuming (3 > and a > for dcfinitcness, the function if (y) must van- 
ish both as y — > —cyo — iO and as y — > — oo, where <£(y) diverges. Setting 
t = —l/(y + cyo), the integrals involved in the above condition can again be 
performed in terms of modified Bcssel functions. Our end result reads 

Q = D ww G{x), (5.8) 

where the scaling function 

K' (x) 

C(x)=-*4i (5.9) 

is identical to that of the supersymmetric monolithic disorder (see (4.24)); we 
will give an explanation of this observation at the end of Sec. 9. Its argument 
reads x = y/&fi/D ww , i.e., 

_ y/a[u-a + 2(D WW - w)cy - , ■. 

x — — . [O.W) 

Uww 

We recall that c = ±1, while yo and v have been defined in (5.1) and (5.6). 
5.2 The partly correlated case (c 2 < 1) 

Let us now turn to the generic partly correlated case where c 2 < 1, so that 
b > 0. The polynomial 

Q(y) =4D ww (y 2 +2cy y + y 2 ) (5.11) 
has two simple zeros at the complex conjugate points 

3/i = (-c - ib)y , ?/2 = {-c + ib)y , (5.12) 
and a double zero at infinity. Equation (1.49) yields 

H(y) = (y- yir-\y - y^- 1 e* (v \ (5.13) 



where 



and 



b + ic i(u — a) + 2w(b — ic)y a + ia(b — ic) 2 y 2 
Vl = ~2b~ + 4D ww by Q ' 

= b-ic i(q -u) + 2w(b + ic)y - ia(b + \c) 2 yl 
V2 2b 4D ww by a ( ' ' 



3.3 



Equation (1.51) reads 

2D ww K\y) = (ay-w + Q)(y - yi )~^ (y - y 2 )-* e"*^. (5.16) 

Assuming Rei^i < 1 for definiteness, the function K{y) must vanish both as 
y — > yi — iO and as y — » — oo. The change of variable 

t=-^-, (5.17) 
2/1 ~ 2/2 

mapping j/i onto and yi onto — 1, leads us to the condition 



Jo 



{ [D ww x(2t + 1) + D ww (y x + u 2 - 1) - Q] 

xt- /1 (t + l)- V2 e- xt }dt = 0, (5.18) 

with 

_ iabyp _ ia^/ D ww D U u - D wu , . 

x - D 2D 2 ■ ^ Ay > 

'-'WW A1J WW 

The integrals involved in the above condition can be expressed in terms of the 
Whittaker function [44, vol. I, (18), p. 274], [47, (9.220), p. 1059] 

Wi, m {x) = — n / t m - 1 —* (t + l) m + l +2 e ~ xt dt (5.20) 

r(m - I + 5 ) Jo 

and of its derivative with respect to x. The end result is 

Q = D ww G(x), (5.21) 

where the argument x has been defined in (5.19), and the scaling function 
G(x) reads 



G(x) = l-2x-^-. 



(5.22) 



In this expression, 



i = V - -2) = 4 + i M ' Q ' 2WC ;° + a(1 " 2c2)?/o2 , (5.23) 
2 V ; 26 4£> t0t0 6i to 

1 ci \ W+acyo 1 /_ aD wu \ 

m=^l-ui-U2) = -— = - l w+ \ (5.24) 

so that Z is imaginary and to real. The numbers c, 6 and yo were defined 
in (5.1). 

In Sec. 9, we will establish the relation between the result (5.22) and a 
result of [48] for a continuous model: the Schrodinger equation with mixed 
disorder. 

Interestingly, a characteristic function with a similar structure — although 
not exactly the same — appears when considering a product of elements in the 
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two-parameter subgroup of SL(2, M) obtained by setting w — in the Iwasawa 
decomposition: 

M=( C ° Sa - S[na ) (I "V (5.25) 
y sin a cos a J \ 1 J 

When a and u are two independent exponential variables, the matrices M n 
are not close to the identity matrix, and yet it is possible to obtain the in- 
variant measure / and the characteristic function fl in closed analytical form. 
Such a product of matrices corresponds to a quantum-mechanical model for 
(^-impurities such that the spacing between consecutive impurities and the 
strength of each impurity are independent and exponentially distributed (see 
Sec. 9.1). This model was first solved by Nieuwenhuizen [10] and studied again 
in a broader context in [21]. 

The differential equation satisfied by the Whittakcr function [47, (9.222), 
p. 1060] 

4x 2 W" m {x) - (x 2 - Alx + 4m 2 - l)W Lm (x) = (5.26) 
translates into the following Riccati equation: 

G(xf - 2xG'{x) = x 2 - Alx + 4m 2 = ^ + ^ Dwu (5.27) 



WW 



for the scaling function G(x). Setting 



A{x) = \J x 1 - 4/x + 4m2 = vV + ^A^ (g 2g) 

Dww 

and solving (5.27) iteratively yields the large- A{x) asymptotic expansion 
, , x(x-2l) x \x 3 + A(l 2 - 4m 2 )x + 16lm 2 ] 

G ^- A ^ +A Mxr- A — + -- (5 - 29) 

When this is expressed in terms of the covariances, we recover the weak- 
disorder expansion results (3.15) and (3.19). 



6 Independent disorder with zero mean 

In this section we deal with the case where the three random variables a n , 
w n and u n are independent and have zero mean. The non-zero parameters 
are therefore the three variances D aa , D ww and D uu . In this situation the 
matrices M n fluctuate around the unit matrix according to 

M n -I^( Wn ""-""V ( 6 .i) 
\a n -w n J 

The mean square norm of the difference between M n and the unit matrix reads 
e 2 = E(|M„ - if) = 2D aa + 2D WW + D uu . (6.2) 
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The polynomial Q(y) reads 

Q(y) = D aa y 4 + 2{D aa + 2D ww )y 2 + D aa + D uu . (6.3) 

Let us henceforth consider the generic situation where none of the variances 
vanishes. This is the first instance where Q{y) has four distinct zeros. Further- 
more we have the polynomial identities 

a(z)a'{z) = 2v(z) and Q'(y) = 2R(y). (6.4) 

The first of the above identities ensures that the condition (2.20) is auto- 
matically fulfilled. The stationary current j therefore vanishes, and the char- 
acteristic exponent fi is real. 

The second of the identities (6.4) yields 

Equation (1.51) then reads 

The next step consists of investigating the zeros of Q(y). First, Q(y) is an 
even polynomial. Setting t = y 2 , we have 

Q(y) = D aa t 2 + 2(D aa + 2D ww )t + D aa + D uu . (6.7) 

The two zeros of this expression read 



D aa + 2D WW + y/(D aa + 2D WW ) 2 - A 2 
1 = D ' 



*2 = pj , (6.8) 



with 



A = VD aa (D aa + D uu ). (6.9) 
Two regimes are to be considered separately. 

Regime I (A < D aa + 2D WW , i.e., D uu < 4D ww (D aa + D ww )/D aa ) 

In this first regime, t\ and t% are real and obey t\ < t 2 < 0. The four zeros of 
Q(y) are ordered as follows along the imaginary axis: 

2/1 = -iV-h, 2/2 = ~W-h, 2/3 = iv/-*2, 2/4 = i\Z^h- (6.10) 

The function K{y) must vanish at both zeros of Q{y) which lie in the lower 
half-plane, i.e., y\ and y 2 , where H{y) diverges. We thus obtain the condition 

"»^ + D + 2fl dy = a 



Jyi 
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y = — - — , - — - — in — 



An integration by parts, corresponding to the choice 

g(y) = — — (6.12) 
y - 2/4 

of the gauge function (see (7.11)), removes the term in y 2 in the numerator of 
the integrand. The change of variable 

yi { m - m) + 2/4(2/2 - 2/1) sin 2 e 

2/4 - 2/2 + (2/2 - 2/i) sin 2 6 
allows one to express the remaining integrals in terms of the Lcgcndrc complete 
elliptic integrals 

/•t/ 2 AO fx/ 2 , 

K(k)= - E(fc)=/ Vl-fc 2 sin 2 0d6», (6.14) 

Jo VI - k 2 sin 2 6 Jo 

where the square elliptic modulus is the cross ratio 

,2 l \ (2/1 - 2/2) (2/4 - 2/3) 

k = 2/1, 2/4; 2/2, 2/3 ) := 7 T7 r- 6.15 

(2/1 -2/3) (2/4 -2/2) 

The zeros y\ , . . . , j/4 can be eliminated from the resulting expression by 
means of (6.8) and (6.10). Some algebra leads us to the following end result 

O = \ ({D aa + 2D WW + - (D aa + A)) , (6.16) 

where the elliptic modulus reads 



7 &aa ~\~ ^D ww Z\ 

k - y D aa + 2D WW + A> [bA{) 

with A being defined in (6.9). 

Finally, the differential identity [44, vol. II, p. 322], [47, (8.123.2), p. 907] 

E(jfe) = (l-k 2 )(K(k) + kK'(k)) (6.18) 

and the linear differential equation [47, (8.124.1), p. 907] 

k(k 2 - l)K"(fc) + (3fc 2 - l)K'(fc) + fcK(Jfe) = (6.19) 

allow one to recast the above result as 

Q = AG(k) + D ww , (6.20) 

where the scaling function 



G(k)^k K '^ 



K(k) 1 - k 2 



(6.21) 



obeys the Riccati equation 



G{kf + kff(k) _ _ D^-^D^, D„) (6 22) 
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Regime II (A > D aa + 2D WW , i.e., D uu > 4D ww (D aa + D ww )/D aa ) 

In this second regime, t\ and t 2 form a complex conjugate pair. The four zeros 
of Q(y) are still given by (6.10). They now sit at the vertices of a rectangle 
in the complex plane, with y\ and j/2 being in the lower half-plane, and their 
respective opposites 2/4 and 7/3 in the upper half-plane. The elliptic modulus k 
is accordingly found to be purely imaginary. The characteristic exponent ft 
can be obtained as the analytical continuation of the result (6.16) to imaginary 
values of k. It is however advantageous to perform the change of modulus [44, 
vol. II, Ch. 13] from k 2 to 

12 k2 ( x (2/2-2/1X2/4-2/3) , ROQ x 

k =-z 72 = (2/2, 2/4; 2/1, 2/3) := 7 77 7- 6.23 

(2/2 - 2/3)(2/4 - 2/i) 

Our result (6.16) thus becomes 

n = A^--\{D aa + A), (6.24) 
K(fc) 1 

with 



k =\l 2A ' (6 - 25) 

The above result can be recast as 

Q = AG{k) + D WW , (6.26) 

where the scaling function G(k) = G(k), i.e., 



G{k) = k(l - k z ) 



2l K'(fc) 



K(fc) 



(6.27) 



obeys the Riccati equation 



G(fc) 2 +fc(l-fc 2 )G'(fc) = k 2 {l-k 2 ) = DaaD ™ {D + D WW ) (g 28) 

4-Daa (-Daa ~4~ Uuu ) 

The rightmost side of (6.22) and (6.28) is the same rational expression of the 
variances, which is negative in the first case, and positive in the second. 

In the borderline situation where D uu — 4:D ww (D aa + D ww )/D aa , so that 
A = D aa + 2D WW and k = k = 0, all the above expressions consistently yield 

n = D ww . (6.29) 



The above general results (6.16), (6.24) can be made more explicit in the 
special cases where one of the variances vanishes. 
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6.1 Supersymmetric and scalar disorder (D aa = 0) 
In this case, (6.16) leads to the simple result 

Q = 0. (6.30) 

Indeed D aa = yields A = and k = 1; the result therefore follows from the 
fact that E(l) = 1, whilst K(fc) diverges to infinity as 

K (*> Mln 7r=P (6 - 31) 

in the fe -> 1 limit [47, (8.113.3), p. 905]. 

6.2 Distance and supersymmetric disorder (£)„„ = 0) 

In this case, we have A — D aa , and so the results (6.16), (6.17) become 



f2 = (D aa +D ww )^§--D aa and k = J - D ™ w . (6.32) 
K(fc) V D aa + D ww 



Further simplifications arise in various regimes. 

For D ww <C D aa , i.e., k — > 0, the expansions [47, (8.113.3), p. 905] 

K(*) = !5>* a » E(fc) = ^a„ T ^-, a„ = |gi (6.33) 

n>0 n>0 1 ' 7 

of the elliptic integrals yield 



n D 2 D 3 

2 16£> aQ 32£)2, 



In the opposite regime (-D^u, 3> D aa ), the behaviour (6.31) yields 

2D 

Q _^ A-LSww /g or\ 

ln(16^/D 

This inverse logarithmic singularity is reminiscent of (4.41). 

When both variances are equal (D ww — D aa ), the modulus reads k = 
l/y/2. In this case, the elliptic integrals can be expressed in terms of Gauss's 
lemniscate constant 

r(i/4) 2 

G = ] 2 = 0.834626 . . . (6.36) 



We thus end up with 



Q = = 0.456946 . . . D aa . (6.37) 
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6.3 Distance and scalar disorder (D ww = 0) 
In this case, the results (6.24), (6.25) read 



n=A W)- 1 * {Daa+A) ' k= y^^- (6 - 38) 

Further simplifications occur in certain limits. 

For D uu <C D aa , the series (6.33) for the elliptic integrals yield 

D qn 2 39 D 3 

8 256D + 2048£>* + ^ 



In the opposite regime (.D uu S> D aa ), the elliptic moduli go to fc = i and 
fc = 1/V2, so that we get 



Q ~ 2^G 2 "" = °- 228473 ' • ' ^D aa D uu . (6.40) 



6.4 Maximal Lyapunov exponent 

The complex characteristic exponent ft has been found to vanish whenever 
two of the variances vanish. This property is due to the fact that the matri- 
ces M n have a trivial composition law in these circumstances: they commute 
among themselves. The above property can thus be viewed as a special case 
of (4.14), (4.43), (4.70). 

It is therefore of interest to look at the maximum of the characteristic 
exponent, for a fixed value of the mean square norm e 2 (see (6.2)). A numerical 
investigation of the result (6.16) shows that fi reaches its maximum 

12 = 0. 170787995... e 2 (6.41) 

for the following values of the variances: D aa = 0.042658 ... e 2 , D ww = 
0.416226... e 2 , D uu = 0.082228. . .e 2 . The corresponding elliptic modulus 
reads k = 0.919798 . . . 



7 The general case 



We now turn to the general case, where the nine parameters of the problem, 
namely the three mean variables and the six covariances, take generic values. 
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7.1 Preliminaries 

The polynomial Q(y) generically has four distinct zeros j/j (i = 1, ...,4). It 
can therefore be factored as 

4 

Q(v) = D*a]l(v-Vi)- (7.1) 
It will also prove useful to introduce the partial fraction expansion 

4 



R(y) 



z ' 7/ — ?/„• 



<2(y) ^ y - y-. 

where the exponents 

a - R (Vi) (73) 

are, in general, complex numbers. Expanding (7.2) around oo to order 1/y, we 
obtain the sum rule 



I> = 2 - ( 7 - 4 ) 

i=l 

Equation (1.49) yields 

4 

H(y) = II(y-wr 04 . (7-5) 
i=i 

so that (1.51) reads 

4 

Daa«"(y) = + 2/2) IJG, - yi ) 04_1 . (7-6) 

i=i 

The polynomial Q(y) = v 2 (y) is manifestly positive for all real values 
of y. As a consequence, no zero of Q(y) can cross the real line. Under generic 
circumstances, two zeros of Q(y), say y\ and j/ 2 : are in the lower half-plane, 
whereas the other two are in the upper half-plane. We choose to label those 
zeros so as to have 

2/3 = 3/2, 2/4 = 2/1, az = a* 2 , a A = a\. (7.7) 



7.2 Calculation of the characteristic exponent 

Let us consider for definiteness the case where the exponents a\ and ai have 
positive real parts. In this situation, the function K{y) must vanish as y = yi 
and y — y 2 , where H{y) diverges. We thus obtain the condition 

/vi JL 
(S(y)+2f2)Y[(y-y^- 1 dy = 0. (7.8) 
»=i 
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The above expression has a form of gauge invariance. This property has 
already been mentioned in two special cases in Sec. 4.3 and 6. It is now time 
to explain it in full generality. Consider an arbitrary regular function g{y). 
Equation (1.49) implies 

d g(y) = Q(y)g'(y) + R(y)g(y) (7 q . 

dy H(y) Q(y)H(y) ' 1 ' ] 

and so 



/ 



-Q(y)g'(y) + R(y)g(y) _ 

Q(y)H(y) v ' 



As a consequence, the expression (7.8) for the characteristic exponent fl is left 
invariant if S{y) is replaced by 

(y) = S(y) + Q(y)g'(y) + R(y)g(y). (7.11) 

The integrations by parts performed so far to simplify the expression for Q in 
two special cases in Sec. 4.3 and 6 correspond to simple rational choices for 
the gauge function g(y) (see (4.52), (6.12)). 

In the general case under consideration, it is advantageous to consider the 
gauge function 

M = - j^, (7-12) 
where y 2 is the endpoint of the integral (7.8). We then have 

S^(y) = ^-+B(y 2 ), (7.13) 

y-yi 

where 

A(y 2 )=Q'(y 2 )-R(y 2 ), B(y 2 ) = \q" \y 2 ) - R> \y 2 ) + S(y 2 ) (7.14) 

are polynomials in y 2 of respective degrees 3 and 2. The condition (7.8) there- 
fore becomes 



rV2 
J VI 



f^A + B (y 2 ) + 2n) f[(y y^ 1 dy = 0. (7.15) 



Now, let us perform the rational change of variable 

t = y_yi^ / 716 x 
yi -y 

which maps yi to 0, y 2 to oo, oo to —1, and j/4 and y 3 to — 1/k and — 1/A 
respectively, where 

yi - yi , Ws - IJ2 ,„ , 

k= , A= . ('-17) 

2/4 - yi 2/3 - yi 
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The integrals involved in (7.15) are thus reduced to integrals of the form [43, 
vol. I, (35), p. 312] 



r(c)r(2 -a-b-c) 
K c r(2-a-b) 
x 2 Fi(l-6,c;2-a-6;a;), (7.18) 



where the argument 



x=l = 2/1, 2/45 2/2, 2/3 := 7 77 r (7.19) 

« (2/1 -2/3 ) (2/4-2/2) 

is the same cross ratio of the four zeros as the square elliptic modulus k 2 
of (6.15), while 



r(a)r(p) 



E 



n>0 



r(n + a )r{n + p) x n 



is Gauss's hypergeometric function. 
We thus get the following expression 



2 



D aa (2/1 - 2/4)(3/3 - 2/2) ai 



2 Fi(l - a 3 ,ai + l;ai + a 2 ;x) 
2 Fi(l - a 3 ,ai;ai + a 2 ;cc) 



+ D aa (2/2 - 2/3)(2/2 - 2/4)(l - 02) - -8(2/2) 



(7.21) 



Finally, the differential relation [44, vol. I, (23), p. 102] 

/3 2 Fi (a, /3 + 1; 7; s) = /3 2 Fi (a, /3; 7; x) + a: ^'(a, /3; 7; x) (7.22) 
allows one to bring the above result to the compact scaling form 



n = I (iWyi - - + - 5(2/1) ) . 
* \ 2/1-2/3 



(7.23) 



where x is the cross ratio introduced in (7.19), and the scaling function G(x) 
reads 



G(x) = x 



2 F/(1 - a 3 ,ai;ai + o 2 ;x) 
2 Fi(l - a 3 ,ai;ai + a 2 ;x) 



(7.24) 



Our end result (7.23) involves, besides basic parameters, the complex ze- 
ros yi of the polynomial Q(y) and the associated complex exponents dj. The 
labeling (7.7) of the complex zeros allows us to recast the argument x as 



x = 



2/2 - 2/i 



(7.25) 



2/2 - 2/i 

This quantity manifestly obeys < x < 1. This condition, together with 

Re(ai+a 2 ) = l, (7.26) 
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ensures that the hypergeometric series entering (7.24) are convergent. The 
above result therefore provides a meaningful expression of the complex char- 
acteristic exponent for generic parameter values. 

The differential equation satisfied by Gauss's hypergeometric function 

x(l - x) 2 F 1 "{a, /3; 7; x) + (7 - (a + /3 + l)x) 2 F 1 '(a, /3; 7; x) 

- a/3 2 F 1 (a,/3;7;x) =0 (7.27) 

translates into the following Riccati equation 

(1 - x)G(x) 2 + x{\ - x)G'(x) 

+ [(a 3 - a\ - l)x + ai + a 2 - 1] G(x) 

= a 1 (l-a 3 )x (7.28) 

for the scaling function G(x). Equation (7.28) provides a useful check of the 
weak-disorder regime, where all the exponents a, become simultaneously large. 
To leading order, we obtain 

G(x) = 2(ayiy 3 - W(t/i + y 3 ) - u + a ± n(yx - y 3 )) 
D aa {yi-y 3 ){yi-yA){y2-y 3 ) 

Choosing the lower sign, we recover after some algebra the expected result 
[} = fi. 

It is worth emphasizing that both the sum rule (7.4) and the choice (7.12) 
of the gauge function are essential for the above reduction to be worked out. 
For instance, the integrals involved in the original expression (7.8), with S(y) 
being a quadratic polynomial, would involve an extra multiple pole at t = — 1, 
and could therefore not be reduced to the form (7.18). Other admissible choices 
of the gauge function g(y) and/or of the rational change of variable from y 
to t correspond to permuting the complex zeros y^. The effect of these discrete 
transformations is to change the argument x into one of the six possible values 
of the cross ratio of the four zeros: 

f 1 1 x x — 1 1 

\x, 1-x, -, , -, \. 7.30 

L XL — XX— 1 X J 

Only the first two possibilities ensure that the argument is in the range < 
x < 1. 



7.3 The case of a vanishing stationary current 

To close, we return to the case where the stationary probability current (inte- 
grated density of states, rotation number) j vanishes, so that fl is real. 

The condition for this to happen has been derived in Sec. 2.2. Using (1.45), 
Equation (2.20) can be recast as 



m Az=a , (7 . 31) 
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Finally, the partial fraction expansion (7.2) allows one to perform the above 
integral explicitly, along the lines of (A. 5). The result is i7r(a 3 + a 4 — a\ — a 2 ). 
Using the relation (7.4), we are thus left with the condition 

01+02 = 1. (7.32) 

In other words, as a consequence of (7.4) and (7.7), the exponents a\ and a 2 
obey in general 

Re( ai +a 2 ) = l, (7.33) 

while the condition (2.20) for the stationary current j to vanish is equivalent 
to the extra condition 

Im (en + a 2 ) = 0. (7.34) 
The resulting expression for the scaling function, 



G{x)=x ggfa4lM , 



(7.35) 



is then manifestly real. It can be checked that the full expression (7.23) for 
the characteristic exponent is also real in this case. 

The condition (7.32) can be expanded in terms of the three mean Iwasawa 
parameters and of the six covariances by eliminating the zeros yi. We thus 
obtain a formidable homogeneous polynomial equation of degree six which 
consists of 335 terms. 

If we restrict the analysis to the situation where the mean parameters a, 
w and u vanish, the condition (7.32) reduces to a homogeneous polynomial 
equation in the covariances with only 22 terms. Its explicit form 

39 D 4 D 2 + 16D D 2 D D 2 — 32 D D 3 D D 

3 on !n ion \rfi n 2 



wu 



^D aa D au D wu 4- SD ww D uu D aw D wu 8D aa D ww D au D aw D wu 



' -^-^uu-^aw ' OLJ ww ctu 1 ^ olw aa wu 

I 2 



~\~ (^D ww -\- AD aa D ww D aOL D uu ) 
x (D aa D 2 wu + 2D aa D aw D wu - D uu D 2 aw )D au = (7.36) 

is however not very illuminating. 

To close, let us mention that the above condition does not directly apply 
to the particular cases studied so far, as they correspond to multiple zeros, 
whereas it was implicitly assumed in the derivation of (7.32) that the four 
zeros of Q(y) were simple. 

For independent disorder with zero mean (Sec. 6), we have shown that j 
vanishes identically This case is very special, as the polynomial identities (6.4) 
imply ^ = 1/2 for i = 1, . . . , 4. 



45 



8 A limiting case: Hyperbolic Brownian motion 

In the previous sections we have shown how to compute the characteristic 
exponent of an infinite product of random matrices in the continuum regime 
where the matrices are close to the identity. This approach provides a rather 
complete understanding of one-dimensional soluble random potentials. The 
knowledge of the Lyapunov exponent brings valuable information on the de- 
gree of localization of the wave functions of an infinite disordered sample. 
Another observable, probably more suited for finite samples, is the reflexion 
phase acquired by a particle which is scattered on the boundary of the sample. 
In the limit of a semi-infinite system, this phase converges to a random variable 
with a non-trivial distribution which has been studied by several authors [41]. 
In particular Barnes and Luck [12] have shown how to relate it to the invariant 
measure of the Riccati variable. 

In this section we will return to this problem in the case of general point 
scatterers. This problem can be reformulated in geometric terms as a random 
walk on the hyperbolic plane. For certain choices of parameters of our model, 
this random walk converges in the continuum regime to a hyperbolic Brow- 
nian motion [49]. The purpose of this section is two- fold. On the one hand 
we demonstrate by an explicit calculation that hyperbolic Brownian motion is 
indeed included in our general scheme. Moreover we will use this correspon- 
dence to study the convergence of the Lyapunov exponent to its bulk value as 
a function of the size of the system. 

We consider a disordered sample of finite length L. We assume that the 
potential which consists of n point scatterers has its support on the interval 
[0,L]. In this setting we now consider the following scattering problem 

• For x ^ 0, the wave function tp(x) = c~ lkx +r(k) e ikx represents an incom- 
ing wave which is partially reflected by the sample. 

• For x > L, the wave function is a transmitted wave ip{x) = t(k) e lkx . 

We have 

/V(oA _ n -i(^'(L)\ 

U(o)J _J7 " UwJ' ( } 

where 

i7- 1 :=M 1 - 1 M 2 - 1 -.-M- 1 (8.2) 

and the M n are given by (1.13). Note that the corresponding Riccati variable 
may be expressed in terms of the scattering data 

z(0) V/(0) 1 - r(k) z(L) _ y/(£) _ 

k - ki/>(0) l + r(k) k ~ kiP(L) [ ' 

Set k = 1 for definiteness. The scattering problem by point scatterers can be 
reformulated as an iteration of Moebius maps 

M (z) = (8.4) 

mn - m 2 iz 
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(see (1.6)) acting in the upper half-plane of the complex variable z = x + iy. 
The Riccati variable z(0) can be expressed in terms of the reversed iterates of 
the point z = i 

z(0) = z n {i) := (Mr 1 ° M2 1 o • • • o M^ 1 ) (i). (8.5) 

The backward sequence {z„(i)} can be viewed as the successive approximants 
of a random continued fraction which converges to a random limit for n — » 00. 
Standard results [26] show that it converges to a real random variable x = 
linin^co z n (i) whose density is given by the invariant measure f(x). 

In physical terms, this expresses the fact that, for a system of large size, 
the incoming wave is entirely reflected (r = e I(5 ), implying z(0) = cot(<5/2). 
The random orbit {z„(i)} defines a random walk on the complex plane which 
can be conveniently described in geometric terms using standard tools of hy- 
perbolic geometry which we briefly recall. 



8.1 The Poincare upper half-plane 
Defined as 

H = {z = x + iy I x, y e R,y ^ 0} , (8.6) 
it is globally invariant under the action of 

M: =(T T 2 ) eSL ( 2 ' R )- ( 8J ) 

^ m 2 l 17122 ) 

The hyperbolic distance between two points z, z' G H, given by 

coshd(z, z 1 ) = {X - X>) 1 + f 2 + V ' 2 , (8.8) 

%yy 

is invariant under the group action: 

d(M(z),M(z')) = d(z,z') VAfeSL(2,R). (8.9) 
Using the identity 

2 coshd(i,M(i)) = \M\ 2 (8.10) 

and the group invariance (8.9) one can compute the distance between two 
successive iterates of the random walk, 

coshd(?„(i),z„_i(i)) = coshd(i,.M„(i)) = ^|M„| 2 , (8.11) 

which shows that the step length distribution only depends on the subgroup 
corresponding to 

' e w \ / 1 m\ 

' ! ' ' (8.12) 



c- w J \ 1 
in the Iwasawa decomposition. 
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8.2 Continuum limit 

When the matrices M n are close to the identity, the recursion relation for 
T — TT^ 1 i e 

T n+1 -T n =T n (Af- + \ - /), (8.13) 

turns into a stochastic differential equation on the group SL(2,R) whose gen- 
eral form we proceed to discuss. 

To this end, it is helpful to work with the finite- dimensional representation 
of SL(2,R) which associates to an element M the operator corresponding to 
matrix multiplication by M. Introduce the infinitesimal generators X i} i £ 
{a 7 w,u}, defined implicitly by 

p ax a _ /cos a -sinoA w x w _ f e w \ uXu _ ( 1 u 
\&ma cosa ) ' ^ e~ w J ' ^0 1 

(8.14) 

One can show that (M n } 1 — I) will converge to a Lie algebra-valued white 
noise with drift. Therefore, replacing the index n by the continuous variable t, 
the recurrence (8.13) goes into 

dT 

- d J= Tt E [A<*i + *i %(*)], (8.15) 

where the r]i (t) are three white noises with some correlation matrix and the \ 
are drift terms expressible in terms of the means a, w 7 u. 

In the following we will analyse a particular case and show how it fits into 
our general scheme. We assume that the non-zero parameters of the model are 
D ww = 1/4, D uu — 1, X w = w = —e/2. Therefore 

+ X im {t)) = - £ -X w + ^X wVw (t) + X uVu (t). (8.16) 

i£{a,w,u} 

In our case, the Iwasawa decomposition of T t is of the form 

T.- (J f) f ^ J_ V («■») 



/y(t) , 

The action on z = i yields the following process in H: 

z(t) = T t (i)=x(t)+iy(t). (8.18) 

The Stratonovich differential equation on SL(2,R) (8.15) is then lifted to a 
Stratonovich differential equation on H: 

x(t)=y(t)r] u (t) and y(t) = -ey(t) + y(t) Vw (t). (8.19) 

By integration we obtain 

y (t) = c- £t+B{t \ where B{t) := f dt' r] w (t'). (8.20) 

Jo 
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The process y(t) is therefore simply the exponential of a linear Brownian 
motion with drift. In a suitable clock r(t) the process x(t) is also a Brownian 
motion [50,51]: introducing the time transformation t — > r(t) defined by 



dt 



= y(t) 2 , i.e., T (t) = f dt'e- 2 " t,+2B ^\ (8.21) 
Jo 



and using r? u (r(i)) ( '= } Vu(t) / \f^[p- , we get 

^fy = ^(*))> ( 8 - 22 ) 

which gives 

/■r(t) 

x(t) = / dr'^r'). (8.23) 
Jo 



Using these results we can compute the Lyapunov exponent as well as the 
invariant measure f(x) and show that they do agree with the results obtained 
from the general formulae (1.32) and (1.39). In fact, since we know explicitly 
the joint law {x(t), y{t)} 7 we can even discuss the case of a finite system and 
obtain the convergence to the semi-infinite system in terms of a central limit 
theorem. Using (8.17) gives 

lniT t| - £t /2 (lg) l 

where N(0, 1) is the standard normal variable. The Lyapunov exponent is there- 
fore 7 = e/2. 

The invariant measure may be derived from the limiting law of (8.23). 
The exponential functional of Brownian motion (8.21) converges to a random 
variable on (0, oo) with probability density 

I \ l+e 



Therefore 



One can check that this is indeed a zero-current solution of (1.32) for the 
choice of parameters D ww = 1/4, D uu = 1, w = —e/2. 

For e = 1/2 the invariant measure is given by the well-known Poisson 
kernel which is known to be the exit law of Brownian motion on H. 
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9 Disordered quantum-mechanical interpretation 

9.1 One-dimensional quantum mechanics with point- like scatterers 

As mentioned in the introduction and discussed at length in [21], any product 
of 2 x 2 matrices can be mapped to a Schrodinger equation of the form 

Hip{x) = Eip(x) with E = k 2 , (9.1) 

where the Hamiltonian involves a potential combining two random functions 

H = —^+W 1 (x) 2 -W{(x) + W 2 (x). (9.2) 
The free Hamiltonian 

H ° = '^ ™ 

is related to the compact component of the Iwasawa decomposition, i.e., to 
the rotation matrix 



cos a —sin a 
sin a cos a 



(9.4) 



On the other hand, W\ and Wi are related to the Abelian and nilpotent 
components, i.e., to the diagonal and upper triangular matrices 

6 c l) and (J fj, (9.5) 

respectively. As Equation (1.12) makes clear, the product (1.2) plays the role 
of a transfer matrix acting on the wave function and its derivative. The 
precise correspondence between the product of matrices and the quantum- 
mechanical model is obtained by letting t n :— a n /k be the distance between 
two consecutive scatterers, so that the position of the nth scatterer is given 
by x n := J2p=i tp- (Of course, this only makes sense if a n > 0.) Set 



WAx) = lim w n 5(x — x n + e), (9.6) 

n 

W 2 (x) = lim y^v n 6(x - x n - e), (9.7) 

n 

where v n = ku n and e -4- + . Then the Hamiltonian (9.2) describes an irregular 
lattice of double impurities, in the terminology of [21]. The case where the 
impurity strength u n is distributed according to an exponential law and w n = 
was considered in [21,10], while the converse problem where u n = and w n is 
exponentially distributed was solved in [21]. 



50 



9.2 Continuum limit for D aa = 

In this paper, we have considered the limit of small random parameters a n , 
u n and w n , retaining only the statistical information contained in the first two 
cumulants. Within this scheme, the constraint a n ^ can only be satisfied 
by setting D aa — 0, so that the distance l n = I = a/k between consecutive 
impurities is the same for all n, and the model describes a regular lattice of 
point scatterers. 

The continuum limit corresponds to a high density of weak scatterers. In 
this limit, (9.6) and (9.7) approximate two correlated Gaussian white noises 
such that 

E(Wi(a;)) = w/£ = kw/a, 
E(W 2 (x)) = v/£ = k 2 u/a, 

Cov (Wi(x) Wj(x')) = j dj 8(x - x') = | dj S(x - x'), (9.8) 

where we have introduced the notation Cov(xy) = E(xy) — E(x)E(y). The 
non-negative covariance matrix is 

/ D ww D wv \ / Dww k Dwu \ /i~v (\\ 

Id n i ko k 2 n r ' 

It is useful for future reference to write down explicitly (1.32) (or (A.l)) in the 
case D aa = 0: 

1 



^D ww z -\- 2,D WU z -\- D 



+ [a(l + z 2 ) + 2(D ww -w)z-u] f(z)=j. (9.10) 



2 — uu 



9.3 Monolithic scalar disorder (Sec. 4.1) 

We give the quantum-mechanical interpretation of the results derived in the 
case where D aa = D ww = 0. It is convenient to introduce the parameter 

a 2 := D vv /t = k 3 D uu /a, (9.11) 

in terms of which the two parameters of (1.60) read 

/3 = £(a 2 /2) 1/3 and n = l^{w/t) 2 + v/£ - k 2 . (9.12) 

We recover the result of Halperin [33,31]: 

where the scaling function G{x) is given by (1.64). 
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9.4 Monolithic supersymmetric disorder (Sec. 4.2) 

For the case D aa = = D uu = 0, we introduce the new parameter g := D ww /£ 
in terms of which the parameters in (4.25) may be expressed as v — K(Wi)/g = 
w/D ww and x = ^^Jv/l — k 2 . Then 

i a _, a (jmEzy (9 . 14) 

where G{x) is now given by (4.24). We have thus recovered the result of [45]. 



9.5 Monolithic distance disorder (Sec. 4.3) 

In this case, the random variable a is not necessarily positive. We may never- 
theless give a quantum- mechanical interpretation in terms of random point-like 
scatterers by expressing the phases a n as a sum of two parts: a deterministic 
positive part a describing the evolution of the wave function in the interval 
between two consecutive impurities separated by a distance i, and a random 
part 5a n modeling the effect of some point-like scatterer. The physical inter- 
pretation of the latter is as follows: an electron incident on such a scatterer 
only experiences forward scattering and receives the phase 5a n . With this in- 
terpretation, only the averaged phase carries an energy dependence: a = k£ 
while D aa is independent of k. The two parameters (4.58) entering the scaling 
function (4.57) take the form 

v — 2k 2 £ 1 I 

X= ^kD~ ^ x =M^^ 2kW 1 2+lj2 - ( 9 - 15 ) 



9.6 Mixed case (Sec. 5) 

The study of Equation (9.2) when the two Gaussian white noises W\ and W2 
are uncorrelated (D wu = 0) has been carried out by Hagcndorf and Texier 
in [48]. The correlated case (D wu ^ 0) was also considered in appendix A of 
that paper. Following the idea proposed there, we introduce a linear combina- 
tion of the two noises: 

<l>{x):=-Wi{x) Dir " 



9 D 

U 



V(x) : =lW 2 {x)- = )-^ (W^x) = ) . (9.16) 
V a J D ww \ a J 

From now on we set E = k 2 = 1 for simplicity. We define the three parame- 
ters u, g and a 2 by 

w D wu 
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D„ 



a 



ff2 = 1 D UU D WW - Dg„ 



We may check that the two Gaussian noises V and are uncorrelated: 



E^(x)V(x'))=0, 



(9.18) 
(9.19) 

(9.20) 



and characterized by 



E(<f>{x))=vg, 
E(V(x))=0, 

Coy{<f){x)(j){x'))=g5{x-x'), 
E(V(x)V(x')) = a 2 S(x-x'). 

By using the identity 

W^xf - W[{x) + W 2 (x) = ^{xf + 4>'{x) + V(x) + E , 

with 

_ u _ wD wu _ / £> ml 
a aD ww \2D W 



(9.21) 

(9.22) 
(9.23) 



we may rewrite the Schrodinger equation (9.2) in terms of the new random 
functions 



dx 2 



+ 0(.t) 2 + ^(x) + V(x) +E =:H + E . 



(9.24) 



The spectral and localization properties of the Hamiltonian H were analysed 
in [48] . The analysis uses the Riccati variable 



(9.25) 



From the Schrodinger equation Hip = Eip, where E := E — Eq, we see that 
z(x) obeys the Stratonovich stochastic differential equation 

-^z(x) = -E - z{xf - 2z{x) 4>(x) + V{x) 



( = } -E - 2gvz(x) - z{xf + ^jAg ~z{xf + a 2 v (x), (9.26) 

where rj(x) is a normalized Gaussian white noise. This equation leads to a 
Fokker-Planck equation for the distribution of the Riccati process z(x). Its 
limit probability density / satisfies the inhomogeneous first-order differential 
equation 



N(E) 



E + 2(v + l)gz + z 



z'2 



2 ^ 2 + ^)di 



(9.27) 
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where the current N(E) represents the integrated density of states per unit 
length of the quantum Hamiltonian H defined in (9.24). This equation is easily 
solved; its solution was analysed in great detail by Hagendorf and Texier [48]. 
The relation with (9.10) can be established by noting that the Riccati vari- 
able z introduced in Sec. 1.3 of the present article can be expressed in the form 



A comparison with (9.25) shows that 



A, 



z = z — 



2D V 



(9.28) 



(9.29) 



With the help of this observation we can check that (9.27) coincides with 
(9.10), up to a constant term D wu shifting the energy. This additive constant 
originates in the definition of a continuum limit for the double impurity model, 
and in the interpretation of that limit in terms of correlated Gaussian white 
noises. More precisely, it is related to the choice of the order of the 5-peaks 
carrying the weights u n and w ni i.e., the sign of e in (9.6), (9.7). Indeed, if we 
consider the e — > 0~ limit in the latter equations, i.e., we reverse the order of 
the last two matrices of the Iwasawa decomposition (1.25), then (1.32) leads 
to (9.10) up to a constant term 2D WU . The result for the continuum model [48] 
is therefore given by the arithmetic mean of the discrete expressions in the 
e — > + and e — > 0~ limits. 

Recall that the integrated density of states per unit length measures the 
density of zeros of the wave function ip(x), that is also the current of the 
Riccati variable through R. Hence 



j = aN(E). 
In [48] , the characteristic function 

Q(E + i0+) = j(E) - inN(E) 



(9.30) 



(9.31) 



associated with the Hamiltonian (9.24) was found in terms of a (Tricomi) 
confluent hypergeometric function: 



H(E + i0+) = -vg + gt, 



1 



U' (^ + fl.tz + l; g) 
J [/(^±i + ^ + l;£) 



where 



-^/^ and fl:=i(l 4 /r- 




(9.32) 



(9.33) 



(Note that Equation (102) of [48] contains a misprint, namely a superfluous 
v + 1 in the denominator.) In terms of the Whittaker function, we have 



U(a, b; z) = z- fe /V/ 2 W_ a+b/Ub _ 1)/2 {z), 



(9.34) 
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and hence 

W^-^l-*^). (9.35) 

Equation (9.35) has precisely the same structure as (5.22). With the help 
of (9.17)-(9.19), we easily identify £ with (5.19), up to a complex conjugation, 
and v/2 with (5.24). The parameter 9, where E = E — Eq = 1 — Eq, coincides 
with (5.23), up to a shift of E by D wu . 

Finally, let us point out that, interestingly, the limit of correlated noises 
for (9.2) with D wu — ±y/D ww D uu maps to a quantum Hamiltonian (9.24) 
with V(x) — 0; this is the supcrsymmctric model studied earlier by Bouchaud 
et al. in [45]. This explains why (4.24) and (5.9) coincide. Note, however, that 
the limit a 2 — > is highly singular, as one would expect from the expres- 
sions obtained for the characteristic function, (9.32) and (9.35). This limit 
was carefully analysed in [48] . 



10 Summary 

In this work we have carried out a systematic study of the complex character- 
istic exponent fl of products of arbitrary random matrices close to the identity 
in the group SL(2,R). Our main result is that the characteristic exponent ad- 
mits a scaling form in the continuum limit where the three mean Iwasawa 
parameters a, W, u and the six elements of their covariance matrix are small 
and comparable. This result is fully general; it provides a unified approach for 
many results found in earlier works for specific cases. 

In the above continuum scaling regime, the complex characteristic expo- 
nent fl is a homogeneous function of its nine variables with degree one. Inter- 
estingly, the corresponding scaling function is (essentially) given by the loga- 
rithmic derivative of a special function of mathematical analysis. The general 
case involves Gauss's hypergeometric function, while particular cases involve 
other special functions (Airy, Bessel, Whittaker, elliptic). All these special 
functions (and many others) can be obtained from the confluent hyperge- 
ometric one by limiting procedures [40, Table 13.6, p. 509]. In the present 
framework, however, the type of special function pertaining to each case is en- 
tirely dictated by the singularities of the integrating factor H(y) at the zeros 
of the polynomial Q(y) — the complexified diffusion coefficient u 2 (z). These 
singularities are in turn set by the multiplicities of the latter zeros, according 
to the correspondence given in Table 1.1 at the end of the introduction. 



Acknowledgments 

This work was supported by Triangle de la Physique. 



55 



A Derivation of Equation (1.42) 

In order to derive Equation (1.42), it will be convenient to work with the following equivalent 
form of the integrated Fokker-Planck equation (1.32): 

\ s(z) ' Tz [f(z) ,r2g(z)] + f{z) { M ' s{z) + \ c ' [s(z) x s ' (z)] } = j - {A - 1] 

Let R be a large positive number, destined to tend to infinity. Multiply the above 
equation by l/(y — z) and integrate over z from — R to R: 



2 J-R y - z dz 



2 

dz 



r y ■ 



-/CO 



H ■ g(z) + i c ■ [g(z) x g'( 2 )] j = /"* 2^L. (A.2) 
2 J J-rV-z 



There follows a discussion of each of the terms appearing in this expression. 

Let us deal firstly with the term on the right-hand side. Consider the closed semicircular 
contour C of radius R centered on the origin in the lower half of the complex z-plane. By 
the Residue Theorem, 

<h :27ri, (A.3) 



L 



and so 



We deduce 



ic z-y 

r- R dz r° iRc w de s 

Jr z ~y J-tt Rc w - limy 

f R — _— ► ™3- (A.5) 

J-R y-Z R^oa 

Next, consider the left-hand side. Since the entries of g(z) are polynomials of degree at 
most two, Taylor's Theorem says that 

g(z) = g(y) -(y-z) g'(y) + ^l^L g" (w ). (A.6) 
The first term on the left-hand side of (A.2) may therefore be developed to yield 

2J- R y-z dz 

1 rR Hz ri 1 r R d 

= g(y) ■■= "T [/(*) <T 2 g(z)} -g\y). dz- [/(*) <T*g(z)} 

2 J-R y — z dz 2 J—r dz 

+ ^ ' \ P <MW -*)-£■ [/C0<7 2 g(0] (A.7) 
2 2 J — r dz 

The right-hand side of this expression is the sum of three terms; we deal with each of them 
in turn. For the first term, integration by parts gives 



s(y) ■ \ 



{ i R r R dz 1 
f(z)<r*g(z) - -f(z)<r 2 g(z)\. (A.8) 

[y-z _ H j- R (y-zy j 



We now let R — > oo to obtain 



1 r°° dz 

r^f(z)cr"g(z) 

dz 



/co 



- \ g(y) ■ ^g(y) [°° —^—/(aJ + ig^.aVto) [°° ~ 
2 J-oo (y - 2 J_ ao y - 

- 7 g(y)-<r 2 s"(y) / dz/(z). (A.9) 
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Recalling the definition of the transform F, we conclude that the first term on the right-hand 
side of (A. 7) yields 



2 !/ - z dz 



„— » 5 g(») ■ "ViO] - T 8(W) ■ »V(»)- ( A -i0) 

fl->oo 2 dy 4 

The second term on the right-hand side of (A. 7) can be integrated immediately: 

R 



-g'(y) ■ \ dz A [f (z] a* e{z) ] = - s ' {y) . I f (z) a* e{z) 



-R 



R^oo 



(A. 11) 



by virtue of the Rice formula (2.26). 

The third term on the right-hand side of (A. 7) is the most delicate: 

K ' /fi - 2 l(y) f R dz(y-z)f'(z) 
J-R 



■ <* S( 



s"(y) 



4 



^g'(j/) I" dz(y-z)±{f(z)(y-z)] 
J-R dz 



+ . CT 2 g » (j/) |* d2(j/ _ 2) d_ {j/ _ z)2] (A 12) 

We may use integration by parts for the first and second terms on the right-hand side of 
this expression. We find 

(y) f R dz(y-z)f'(z) — ► I^M. CT 2 g(j/) (A.13) 



4 

and, by using (2.24). 



cr g 



g"(y) 



<r 2 S '(y) f R dz(y-z)±[f(z)(y-z)} 
J-R dz 



■ cr g 



Finally, for the third term on the right-hand side of (A. 12) 
■ ** S "(y) [ R dz(y - z) A [/(*) (j/ - ,) 2 ] 



/oo 
dz(i/ - *)/(*) (A.14) 
-oo 



— o — ■ g 



"(l/) {l/ 2 ^* dz(» -*)/'(«) -2» J\z(y-z)±[zf(z)] 

+ f* dz(,y-z)j- [z 2 f(z)]} (A.15) 



By using integration by parts and (2.24)-(2.25), we obtain 



— o — ' g 



(») J\z(y-z)j- [f(z) (y-zf 



"(y) {?/ 2 - 2y y^°° «/(z)da - 2 A [ z 2 /(z) ] d2 } . (A.16) 
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This completes our discussion of the first term on the left-hand side of (A. 2). 

The second term on that left-hand side is easier to deal with since no derivative of / 
appears and hence there is no need to integrate by parts. By using (A. 6), this term is easily 
developed to yield 

f R R ih s{z) { M ■ s(z) + \ c ' [s(z) x s ' (z)] } 

r^L ' g ^ + \ c ' x g '^] } F ^ ~ M ' g '^ ~ \ c ' x 

+ { \» ■ g"(y) + \ c ■ [s'(y) x g"(y)] } (y - z)f(z) dz. (A.17) 

Equation (1.42) follows after putting all these results together and making use of (1.39) 
for the Lyapunov exponent. 
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